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?  3  ABSTRACT 

A  theory  of  thermoviscoplasticity  is  developed,  which  Is  firmly  founded  on  the 
principles  of  irreversible  thermodynamics  of  internal  variables  developed  previous- 
•  ly  by  the  author.  I  differs  from  previous  theories  In  the  essential  axiom  that 
the  material  senses  oranges  in  deformation  and  temperature  in  a  time  scale  which  is 
an  inherent  material  property.  The  theory  has  thus  been  termed  END0CHR0NIC. 

Material  behavior  is  described  in  terms  of  '’heredity  functions"  akin  to  those  in 
viscoelasticity.  Many  applications  to  the  mechanical  behavior  of  copper  and 
aluminum  are  then  considered,  and  many  facets  of  their  behavior  such  as  yield, 
hysteresis,  Bauschinger  effect,  cross-hardening,  cyclic  hardening,  creep  and  re¬ 
laxation  are  predicted  and  described  analytically  by  the  theory,  with  remarkable 
accuracy. 
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Stcudly  CUktitficolicn 


Final  Report 


This  report  is  a  summary  of  the  work  accomplished  during  the  year  1 970- 1971* 
under  AFOSR  Grant  1916- 70 .  The  work  consists  in  three  main  parts: 

(a)  The  development  of  mathematical  constitutive  theory  for  inelastic 
materials  undergoing  deformation  in  a  non- isothermal  environment. 

(b)  The  application  of  the  theory  to  certain  specific  situations,  where 
experimental  measurements  were  made  by  other  people  in  the  past;  the  purpose 
part  (b)  was  to  establish  the  physical  soundness  of  the  theory  and  to  illus¬ 
trate  its  applicability. 

(c)  The  performance  of  a  set  of  critical  experiments  with  a  view  to  de¬ 
termining  the  material  (heredity)  functions  that  appear  in  the  constitutive 
equation  derived  in  (a).  Knowledge  of  the  material  functions  will  help  in 
their  analytical  representations  necessary  for  the  solution  of  problems 
relevant  to  the  analysis  and  design  of  aerospace  systems. 

The  development  of  the  theory  is  based  on  the  principles  of  Irreversible 
thermodynamics  and  the  axiom  that  the  stress  is  determined  by  the  previous 
history  of  deformation  defined  on  a  time  scale  which  is  itself  a  property  of 
the  material  at  hand. 

The  foundations  are  then  laid  for  a  class  of  constitutive  equations  which 
govern  the  behavio  of  inelastic  materials  when  these  are  subjected  to 
coupled  thcrmomechanlcal  disturbances.  This  is  done  in  the  first  report 
entitled,  "A  Theory  of  Vtscopl  istici ty  Without  a  Yield  Surface  Part  I  * 

General  Theory". 

Applications  of  the  theory  are  dealt  with  In  extenso  in  the  second  report 
entitled,  "A  Theory  of  Viscoplasticity  Without  a  Yield  Surface  Part  II  - 
Appll  cation  to  Mechanical  Behavior  of  Metals". 


2. 


A  number  of  problems  involving  complex  deformat ion  of  metals  under  iso¬ 
thermal  conditions  of  small  strain,  are  given  analytical  solutions  for  the 
first  time.  These  solutions  are  compared  with  experiments  and  such  compari¬ 
sons,  in  all  cares  favorable,  speak  eloquently  for  the  power  of  the  theory. 

In  parallel  with  the  theoretical  effort,  an  experimental  program  was 
undertaken  to  elucidate  a  number  of  points  which  are  likely  to  be  of  import¬ 
ance  in  the  further  development  and  application  of  the  theory.  In  particular 
the  physical  basis  for  the  assumption  of  plastic  incompressibility  was 
tested;  it  was  found  {just  as  Bridgman  did  in  19^8  but  his  work  went  un¬ 
noticed)  that  no  such  basis  exists.  Detailed  discussion  of  this  aspect  of 
the  work  is  found  In  Report  3.  entitled,  "Material  Instabilities  in  the  Ex¬ 
perimental  Study  of  the  Plastic  Compressibility  of  Some  Important  Metals". 

Further,  the  form  of  the  heredity  functions,  that  appear  In  the  theory, 

% 

was  Investigated  experimentally  for  copper  and  aluminum.  Realistic  analytical 
representations  of  these  functions  Is  now  possible  for  the  future  application 
and  further  refinement  of  the  theory.  This  part  of  the  experimental  program 
1$  discussed  in  Report  A  entitled,  "Experimental  Determination  of  the  Heredity 
Functions  of  Copper  and  Aluminum". 

Reports  I  and  2  have  already  been  published  in  the  European  Journal, 
Archivum  Mechanlkt  Sostowanej.  Published  by  the  Polish  Academy  of  Science. 

The  reference  Is:  Arch.  Mech.  Sost.,  (1971). 

In  summary,  visible  progress  (conceptual  as  well  as  analytical)  has  been 
made.  In  this  one  year  tri  the  understanding  of  the  phenomenological  descrip¬ 
tion  of  the  mechanical  behavior  of  metals,  as  a  result  of  this  Grant. 
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ABSTRACT 


Herein,  we  propose  a  mathematical  theory  of  thermo-viscoplasticity 
which  is  a  synthesis  of  experimentally  observed  material  behavior  on  one 
hand,  and  the  concepts  of  irreversible  thermodynamics  on  the  other. 

The  underlying  principle  is  that  the  history  of  deformation  is  de¬ 
fined  in  terms  of  a  "time  scale"  whir*'  is  not  measured  by  a  clock,  but 
is  in  itself  a  property  of  the  material  at  hand. 

The  theory  is  unifying  in  the  sense  that  theories  of  plasticity, 
viscoelasticity  and  elasticity  can  be  obtained  from  it  as  special  cases 
by  imposing  suitable  constraints  on  the  material  parameters  involved; 
furthermore,  it  does  not  make  use  of  the  idea  of  a  yield  surface. 


See tion  1 .  End ochroni c  Theory  of  Viscoplastieity 

In  current  theories  of  plasticity,  to  explain  the  observed  discontin¬ 
uities  in  material  behavior  upon  loading  beyond  the  ’'yield  point”  and  upon 
unloading,  one  has  to  introduce  the  concept  of  a  yield  surface  in  stress 
space  as  well  as  a  "loading  function"  to  distinguish  between  loading  and  un¬ 
loading.  Similarly,  in  the  case  of  viscoplasticity,  the  existence  of  a 
static  stress-strain  relation  and  a  yield  surface  are  assumed  and  the  stress 
increment,  with  respect  to  the  static  value,  is  related  to  the  strain  rate, 
or  more  generally  to  the  strain  history,  by  a  constitutive  equation. 

Howevr,  the  fact  that  the  phenomenon  of  yield  is  usually  a  gradual 
transition  from  a  linear  to  a  non-linear  stress-strain  response,  makes  it 
difficult  to  say  precisely  where  yield  has  occurred,  to  the  extent  that 

♦ 

different  definitions  of  yield  are  used  for  this,  purpose.  Three  such  defin¬ 
itions,  for  instance,  arc  (a)  the  deviation  from  linearity  in  the  relation 
between  soma  measure  of  strain  and  stress  (b)  the  intersection  of  the 
initial  part  of  a  stress  strain  curve  and  the  backward  linear  extrapolation 
of  the  "plastic"  part  of  the  curve  and  (c)  a  value  of  "proof"  stress 
corresponding  to  an  arbitrarily  defined  value  of  "proof"  strain. 

Though,  from  an  engineering  viewpoint,  the  initial  yield  surface  is  not 
overly  influenced  by  the  definition  of  yield,  it  has  been  lound  experimentally 
that  subsequent  yield  surfaces  of  a  strain  hardening  material  are  influenced 
by  the  definition  of  yield  to  an  extraordinary  degree.  (See  Appendix  1).  If 
we  insist  that  the  increment  in  plastic  strain  is  to  be  normal  to  the  yield 
surf'ce,  then,  for  complex  stress  hi  tories,  each  such  definition  will  give 
rise  to  a  different  plastic  strain  history.  Only  one  of  these  can  be  the 


correct  one. 
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The  conceptual  difficulties  that  arc  encountered  by  the  introduction  of 
the  yield  surface  are  completely  circumvented  by  our  theory  of  plasticity 
which  is  developed  on  the  basis  of  the  observation  that  the  state  of  stress 
in  the  neighborhood  of  a  point  in  a  plastic  material  depends  on  the  set  of 
all  previous  states  of  deformation  of  that  neighborhood »  but  it  does  not 
depend  on  the  rapidity  at  which  such  deformation  states  have  succeeded  one 
another*. 

The  independence  of  stress  of  the  rapidity  of  succession  of  deformation 
states  is  achieved  by  introducing  a  time  scale  £  which  is  independent  of  t, 
the  external  time  measured  by  a  clock,  but  which  is  intrinsically  related  to 
the  deformation  of  the  material. 

Of  course  there  arc  many  ways  of  introducing  such  a  time  scale.  However, 
it  appears  almost  mandatory  that  C  should  be  a  monotonically  increasing. func¬ 
tion  of  deformation,  otherwise  two  different  states  of  deformation  could  exist 
"simultaneously"  i.e.  for  the  same  value  of  Furthermore,  a  positive  rate 
of  change  ,  of  the  internal  energy  density  t  with  respect  to  £  could  not  ue 
interpreted  unambiguously  as  a  process  of  increasing  c,  if  dC  could  be  negative. 

A  logical  definition**  for  C  is  then  given  by  the  relation 


*In  the  present  Section  and  in  subsequent  Sections  (with  the  exception  of 
Section  2)  we  shall  assums  that  mechanical  changes  take  place  in  a  constant 
temperature  environment,  such  as  an  isothermal  atmosphere.  The  thermal 
changes  in  the  material  will,  therefore,  be  mechanically  induced  and.  In 
general ,  will  remain  small.  Conversely,  only  thermal  changes  of  this  nature 
will  be  considered  in  this  papsr. 

** Alternative  but  less  general  definitions  have  appeared  in  the  litera¬ 
ture.  For  instance,  Ilyushin^l'  and  later  Rivlin'2'  defined  a  "time"  s  by 

the  relation  ds2  *  dC.dC...  However,  ue  have  found  that  this  definition  is 

O  *3 

too  narrow  to  describe,  quantitatively,  material  behavior  in  the  plastic 
range  as  will  be  discussed  lstsr.  The  effect  of  temperature  on  (,  will  at 
this  time,  be  considered  sufficiently  small  to  be  negligible.  For  s  more 
vague  allusion  to  this  possibility  see  also,  Schapery  (3). 
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d£2 


Pijk£dCijdCk£ 


(3.1) 


where  C...  is  the  Right  Cauchy-Green  tensor  and  is  a  fourth  order  tensor 


which  could  depend  on  C,, ^  .  The  positive  definite  nature  of  d£2  requires 
that  be  positive  definite.  In  the  case  of  small  deformation 


dt2  *  **ijk£dcijdck£ 


where  is  the  small  deformation  strain  tensor  and  could  depend  on 

V 


Actual  materials,  on  the  other  hand,  do,  in  general,  depend  on  the  his¬ 
tory  of  deformation  as  well  as  on  the  rapidity,  or  rate,  at  which  deformation 
states  succeed  one  another.  To  describe  materials  of  this  type  one  may  con¬ 
struct  a  theory  of  viscoplasticity  by  introducing  a  time  scale  C  which  is  re¬ 
lated  to  the  external  time  t. 

It  appears  logical  to  define  C  by  the  relationship 


d C2  *  o2dC2  ♦  B2dt2 


where  e  and  B  are  scalar  material  parameters.  Henceforth  d(  will  be  called 
an  *■  intrinsic  time  measure”,  and  *U),  such  that  jjr  >  0  (()<£<•) ,  will  be 
called  an  "intrinsic  time  scale”. 

In  our  theory,  the  stress  (among  other  properties)  is  necessarily,  a 
functional  of  the  strain  history,  defined  with  respect  to  the  intrinsic  time 
scale,  the  latter  being  a  property  of  the  material  at  hand.  As  a  result  we 
have  called  our  theory  an  endochronic  theory  of  viscoplasticity. 

The  theory  will  now  be  developed  in  a  general  thermodynamic  framework  in 
Section  3.  Before  this  is  done,  however,  the  thermodynamic  foundations  are 


laid  in  Section  2 


4 


Section  7_.  Thermodynamic  Foundations 

The  following  arc  the  fundamental  laws  of  thermodynamics,  which  apply 
to  all  continuous  media  irrespective  of  their  constitutive  properties.  (For 
materials  that  are  solid-like,  in  the  sense  that  they  have  a  memory  of 
their  initial  configuration,  it  is  more  convenient  to  express  these  laws  in 
the  material  coordinate  system  x^).  In  differential  form,  these  are  the 
first  law  of  thermodynamics , 

e  =  (p  /2p)  Tij  C..  -  h\  +  Q,  (2.1) 

o  i]  ,i 

the  rate  of  dissipation  inequality, 

6y  =  (P  /2p )  Ti;5  C.  .  -  *  -  n6  >  0  (2.2) 

013  — 

and  the  heat  conduction  inequality 

-h1*,.  >  0  .  (2.3) 

The  symbols  in  eq.  *s  (2.1),  (2.2)  and  (2.3)  have  the  following  meaning:  Poe 
is  the  internal  energy  per  unit  mass;  pq  and  p  are  the  initial  and  current 
•ass  densities  respectively;  is  the  stress  tensor  in  the  material  coor¬ 
dinate  system  x*;  is  the  right  Cauchy-Green  tensor,  h*  is  the  heat  flux 
vector  per  unit  undeformed  area  in  the  material  system;  p^Q  is  the  heat  supply 
per  unit  mass;  9  is  the  temperature,  y  the  Irreversible  entropy  and  *  and  n 
are  the  free  energy  and  entropy,  respectively,  per  unit  undefomod  volume 
finally  a  subscript  following  a  comma  denotes  differentiation  with  raspect 
to  the  corresponding  material  coordinate.  A  dot  over  a  quantity  denotes 
material  derivative  with  respect  to  time.  To  avoid  repetitious  statements, 
henceforth  we  shall  refer  to  as  the  "deformation". 
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In  the  case  of  dissipative  materials  the  stress,  the  internal  energy 
and  entropy  densities  (and,  therefore,  the  free  energy  density)  of  a  material 
neighborhood  depend  on  the  entire  history  of  deformation  and  temperature  of 
that  neighborhood. 

In  the  theory  of  irreversible  thermodynamics  the  effects  of  history  are 
taken  into  account  by  specifying  that  the  stress  and  free  energy  density  are 
functions  of  the  current  values  of  and  6  as  well  as  n  additional 

independent  variables  q  ,  not  necessarily  observable,  called  "internal  vari~ 

a 

ables".  These  may  be  scalars  or  components  of  vectors  or  tensors  in  the 
material  frame i  whatever  their  geometric  nature  they  must  remain  invariant 
with  translation  and  rotation  of  the  spatial  system  to  satisfy  the  principle 
of  material  indifference.  Thus: 

*  =  *  (C..  ,  0  ,  qa)  (2.4) 

Tij  =  tij  (Cki  ,  e  ,  qa)  (2.5) 

(4) 

It  has  been  shown  elsewhere  that 


s  ?£L  liL_ 

p  3C.. 
o  13 


n  =  - 


ii 

36 


*  *  -  i? 

Furthermore  the  heat  flux  vector  h*  is  a  function  of  6, 


q  i.e., 
a 


hl  *  h1  (6^  ,  «  .  Ck£  ,  q  ) 


(2.6) 

(2.7) 

(2.8) 
6  ,  and 


(2.9) 


subject  to  the  conditions: 


6 


=  0,  hi,.  <_  0.  (2.10  a,  b) 

Finally  eq.  (2.1)  in  conjunction  with  eq.’s  (2.6)  and  (2.7)  yields 


>6,.  « 
*1  =  0 


=  8  (If)*  '  !ja  +  « 


(2.11) 


The  remarkable  property  of  the  above  equations  is  that  they  apply  to  all 
materials  irrespective  of  their  constitution.  This  has  not  been  generally 
recognized.  In  fact  the  constitutive  nature  of  the  material  follows  from 
the  constitutive  properties  of  q  ,  For  example,  in  elastic  materials  o  *  0, 

Q  'O 

whereas  in  viscoelastic  materials  q  are  given  by  a  set  of  differential 

<B 

equations  of  the  type. 


\  s  fa  <qg*  cij»  6)  (2.12). 

The  question  of  how  are  determined  for  viscoplastic  materials  is  consi¬ 
dered  in  the  next  Section. 
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Section  3.  Constitutive  Equations  in  Viscoplasticity 

From  the  right  hand  side  of  eq.  (2.8)  and  the  fact  that  ^  >  0  and 

at 

ds/d(  >  0,  it  follows  that 


-  >_  0  (a  not  summed) 


(3.1) 


where  the  inequality  is  valid  unless  -r—  =0.  It  also  follows  from  inequality 

dqa  „ 

(3.1)  that  -p—  *  q  ,  C. .  and  0  must  be  related  otherwise  r*—  and  >3——  could 
az  a  j]  3q  dz 

0 

be  prescribed  independently  and  in  such  a  fashion,  that  the  inequality  would 
be  violated.  In  this  event  there  must  exist  a  set  of  relations 


3T  *  V'ij  «3-J> 

for  all  0,  where  the  functions  f  are  material  functions. 

o 

It  must  be  noted  that,  as  a  result  of  eq.  (3.2)  qQ  are  indeed  functionals 
of  the  histories  of  deformation  and  temperature  with  respect,  however,  to  the 
intrinsic  time  scale  z  which  is,  itself,  a  material  property. 

Thus,  at  least  formally,  the  constitutive  equations  of  the  endochronic 
theory  of  viscoplasticity  are  now  complete  in  the  sense  that  given  the  mate¬ 
rial  functions  fa,<|i  and  h*  then  for  some  specified  deformation  and  temperature 
histories,  q^  are  found  from  eq.  (3.2)  and  thus  and  n  are  found  from  eq’s. 
(2.6)  and  (2.7)  respectively;  similarly  h^  the  heat  flux  vector  is  determined 
from  eq.  (2.9). 

Ideally,  one  would  like  to  knc  r  the  thermomechanical  three-dimensional 
response  of  a  material  over  the  whole  spectrum  of  mechanical  and  thermal  con¬ 
ditions,  i.e.,  under  all  variations  in  strain,  strain  rate  (or  more  generally, 
history  of  strain)  and  temperature.  However,  such  a  task  would  be  a  momentous, 
if  not  an  impossible,  undertaking;  the  experimental  evaluation  of  the  material 
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functions  involved,  under  wide  conditions  of  strain  and  temperature,  would  be 
impractical . 

Fortunately  the  domain  of  specification  of  design  conditions  is  usually 
limited  in  some  way;  for  instance  usually,  (a)  large  changes  of  temperature, 
fast  rates  of  loading,  but  small  strains  are  prescribed;  or  (b)  small  changes 
in  temperature  and  small  rates  of  loading  but  large  strains  and/or  displace¬ 
ments  prevail.  More  extreme  mechanical  as  well  as  thermal  conditions  are 
rarer. 

It  is  reasonable  to  expect  that  material  behavior  would  be  easier  to 
describe  mathematically  over  a  narrower  domain  of  environmental  conditions, 
where  the  applicability  or  "correctness”  of  such  mathematical  formulation 
would  be  easier  to  check  experimentally. 

In  what  follows  we  shall  consider  situations  in  which  the  strain  id  a 
material  region  R  as  well  as  the  temperature  changes  relative  to  a  uniform 
reference  temperature  6q  are  "small".  To  make  the  above  statement  more  pre¬ 
cise  let  e^j(z')  denote  the  history  of  the  strain  tensor*,  for  zq  <_  z'  <_  z 

where  z  is  some  initial  intrinsic  time.  Set 
o 

» 

||ci;.(z')||  e^iz')}*5  (3.3) 

and  let  the  supremum  of  | |c^j(z')| | be  A. 

Similarly  letdfz')  be  the  history  of  the  temperature  increment  relative 
to  the  reference  temperature  6o  and  let  the  supremum  of  |9(z')|  be  6.  The 
notion  of  smallness  is  made  precise  by  stipulating  that  &«1  and  6<<1. 

Thus,  formally 

6  *  e  ♦  $  (3.4) 

0 


*‘ij  *  *  (eij  -  v 
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i 


n  =  n0  +  x 


(3.5) 


|-5(z')|  =  t 

'  'sup 


(3.6) 


To  complete  the  formalism  let  x  be  the  entropy  change  relative  and  a 
reference  uniform  entropy  no  ,  and  let  denote  the  stress  tensor.  The  re¬ 
ference  state  is  defined  by  the  condition  that  =  0,  <|>  =  0*x  =  0*  qa  *  0. 
Under  these  conditions,  eq.'s  (2.4)  through  (2.9)  and  eq.  (2.10b)  become. 


*  =  *  Ujj  .  5  ,  qa) 

(3.7) 

c  .  it 
-  'll  »«y 

(3.8) 

*  --ft 

(3.9) 

V  =  -  Ua  q*  i  0 

(3.10) 

hi  *  kij  5  *j 

(3.11) 

hi  ^»i  -  0 

(3.12) 

where  is  the  thermal  conductivity  tensor. 

Finally  eq.  (2.11),  becomes 

hi,i  *  e.(H)  Q  ■ 

(3.13) 

It  is  shown  in  Appendix  II  that  |qQ|  and  |$Q|  may  stay  small  in  the  sense 
that  given  two  positive  numbers  6*  and  6g  ,  however  small ,  A  and  <  can  be  cho¬ 
sen  small  enough  such  that  (q^  <  6*  and  | |  <  6®  • 

At  this  stage  one  may  obtain  the  corresponding  equation  for  q0  by 
linearising  eq.  (3.2).  However  it  is  physically  more  meaningful  and,  a  pos¬ 
teriori,  more  rewarding  to  examine  more  closely  the  rate  of  change  of 


irreversible  entropy  y  . 
From  eq.  (3.10), 
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6  -  -  It  >  n 

e  dz  ”  “  3q  dz  -  0 
na 


(3.14) 


It  follows  from  eq.'(s)  (3.14)  and  (3.2)  that  -r-  may  be  expressed  as  a  func- 
dq  az 

tion  of  t  c  and  0  ,  subject  to  the  condition  that  g—  =  0  whenever 

*%  dv 

-  0  for  all  a  .  Thus,  if  we  expand  0^  *  Tavlor  series  and  ignore 

terms  of  order  higher  than  0(6^)*  and  observe  the  inequality  (3.10),  thereby 

eliminating  the  linear  terms  in  the  expansion,  then 


,  il.k 

o  dz  a0  dz  dz 


(3.15) 


Eq.'s  (3.14)  and  (3.15)  are  simultaneously  satisfied  if 


(3.16) 


With  eq.  (3.16)  the  constitutive  description  of  a  viscoplastic  material  is 
now  complete. 

Explicit  Constitutive  Equations 

Explicit  constitutive  equations  for  viscoplastic  materials  under  condi¬ 
tions  of  small  strain  and  small  changes  in  temperature  arc  obtained  by  ex¬ 
panding  ♦  In  eq.  (3.7)  in  Taylor  Series  and  omitting  terms  of  order  higher 
than  0(62)j  *  linear  terms  must  vanish  to  satisfy  the  initial  conditions. 

Before  the  expansion  is  carried  out,  however,  it  appears  desirable  to 
regard  qft  not  as  scalars  but  components  of  second  order  tensors.  This,  as 
will  be  shown,  obviates  certain  difficulties  which  arise  with  the  represen¬ 
tation  of  fourth  order  tensors.  For  instance,  in  Ref.  (4),  we  were  faced 


is  the  largest  of  6*  and  6®.  Also,  6  is  the  largest  of  6,  A  and  A®. 
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with  having  to  assume,  without  proof,  that  a  fourth  order  tensor  C..  .  such 

t  ^  V/ 


that, 


Cijk£  s  Cj ikt  =  Cij£k  =  Ck£ij 


is  given  by  the  series 


„  _  r  aijaak£e 

cijk£  ’  i  >r — 


(3.17) 


where  a.,  are  second  order  symmetric  tensors  and  a  are  scalars.  Problems 

ija  o 

such  as  this  are  obviated  by  giving  the  internal  variables  a  tensorial  char¬ 
acter.  Thus  the  free  energy  density  and  other  thermodynamic  quantities  are 
now  functions  of  ,  9  and  n^  internal  variables  (a  =  1,  2  .  .  .  .  n), 
where  q^  are  symmetric  second  order  tensors,  with  respect  to  the  material 
system  xi  .  In  this  notation,  eq.'s  (3.13)  through  (3.16)  now  read, 


hi,i  eo  39  "  .a  qij  *  Q 

*qij 

,  SL.  .11  lii.  o 

o  dz  .  o  dz  — 

Jqij 


i  *  Plj  m  |a!j  | m 


(3.18) 


(3.19) 


(3.20)* 


3*  .o 

r.  ♦bij«ar“ 


(o  not  sumned) 


(3.21) 


08  *^qk£  at  a  6 

•Expansions  of  the  type  g; —  «nd  Aijk4  qij  reduce  to 


the  above  form.  See  Ref.  5. 
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Furthermore,  in  view  of  my  previous  discussion, 

*  S  ^ijk l  Eij  €k l  *  Bijk£  Cij  qk l  +aCijk£  qij  qk l 

+  D..$e..  +  E?.$q?.  +  *sF$2  (3.22) 

i]  I]  ID  HiD 

Though,  in  principle,  eq. ’s  (3.8),  (3.9),  (3.18),  (3.21)  and  (3.22)  are 
sufficient  for  the  derivation  of  explicit  constitutive  equations,  we  shall 
obtain  these  only  for  isotropic  material,  so  as  to  keep  the  algebra  at  a 
minimum.  For  such  materials 


*ijk*  ’  A1  4lj  6kl  *  *2  4ik  4ji 


EV  *  B1  4ij  V  4  B2  4ik  ‘it 


*  Cl  4ij  \l  *  C2  4ik  4« 


D1J  *  "  4lj 


(3.23  »-f) 


EIl 


kij  *  k  4ij 


*>Ukt  ■  b°  4lj  \l  *  b2  4ik  v 


It  is  worth  noting  that  here  we  consider  materials  which  are  "stable" 

in  the  sense  that  straining  of  the  reference  configuration  under  isothermal 

conditions  will  increase  the  free  energy  density  +  .  Thus  A^^  «nd 

are  positive  definite.  As  a  consequence  ,  C®  and  are  all  positive. 

Ctoitting  superfluous  algebra,  the  coupled  thcrmonechanical  constitutive 

°kk 

equations  take  the  following  form  in  terms  of  the  hydrostatic  stress  o  *  -j-  , 
the  deviatoric  stress  tensor  o ^  ,  the  increment  in  temperature  3  ,  the 


« 
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hydrostatic  strain  ,  the  deviatoric  strain  tensor  e^..  and  the 
entropy  increment  x  ;  in  terms  of  the  above  notation: 

,z  3e 

•ij  =  2  j  y  (z"z,)  3rldz' 

2* 

0 

a  -  j  K( z-z* )  ^dz*  +  j *  D(z-z')  ||,  dz* 

*0  zo 

*X  *  IS  *  f  «-*•>  ♦  f  ««’>  H 

'z_  ;z. 


dz' 


where. 


B®  B®  B®  B?  -p  z 

2  p<z)  =  (A  -  l  )  H(z)  +  l  -?-4  e  ® 


2  -  .a 

®  C2 


a  c: 


B®  B®  B®  B®  -X  z 

K(z),  (A0-S-^)  H(z)  +  X^2.  c® 
C  a  C 


BE  B  E  -X  z 

D(z)  *  (D-  l  44)  H(z)  tI-54  •  * 


o  c: 


a  c: 


-  _  r®  r®  r®  r®  *1.* 

r(l)  .  (r  -  I  £-;£  )  HU)  +  J  .  * 

•  co  •  c. 


(3.24) 


(3.25) 


(3.26) 


(3.27) 


(3.28) 


(3.29) 


(3<a0) 


Ao  ■  1/3  (3Aa  ♦  A2>  ,  Bo»  1/3  (3BX  ♦  B?)  ,  C0  *  1/3  (3^  ♦  Cj)  (3.31) 

(3.32) 


C®  C® 

'*»  '’o 

p  «  -t  ,  x  «  ~ 

o  .a  *o  vo 

b2  bo 


>:  ■  »;  ♦  *>; 


(3.33) 


The  heat  conduction  equation  it  similarly  found  to  bet 


a. k  3.ii  *  tl  1]  '  ft  «*’  *  '.if  f]  «*-*•>  «*• 


where 


C  (z)  =  -  0^  FCz)  (3.35) 

v  o 

and  a  roof  over  a  quantity  implies  differentiation  with  respect  to  z.  The 
lower  limit  z  denotes  the  intrinsic  time  of  the  reference  state. 
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Section  4.  Endochronic  Theory  of  Plasticity  and  its  Relation  to  Present  Theories 

Our  theory  of  plasticity,  which  is  a  rate  independent  endochronic  theory, 
is  obtained  by  replacing  the  time  measure  d(  by  d£.  The  time  scale  now  becomes 
*(0,  but  the  form  of  the  constitutive  equations  remains  unaltered.  In  parti¬ 
cular,  the  "linear"  form  of  our  theory  is  obtained  by  setting 


d*  *  Pijk t  dcij  Atkt 


(4.1) 


where  Pjj)^  is  a  positive  definite  fourth  order  material  tensor.  We  repeat 
the  constitutive  equations  of  the  linear  theory,  in  the  particular  case  when 
the  deforms*...  _«<  is  isothermal  so  that  a  comparison  may  be  made  with  current 
theories.  When  9  £  0,  then  eq.'s  (3.24)  and  (3.25)  become. 


•ij  *  2  j  »<*-*'>  KT5-  da* 
okk  *  3  I  K(s-z')  ~dz' 


<4.2) 


(4.3) 


where  s  «  z(0. 

If  the  material  behaves  elastically  under  pressure  (so-called  plasti¬ 
cally  incompressible)  hen  K(t )  is  a  constant  and  in  this  case 


•toe  ■  5K'kk 


(4.4) 


Whereas, 


■ij  • *  J] » <*-*’  >  d,D(*’  > 


(4.5) 


Let  now  u  (s)  consist  of  a  single  exponential  term  i.e. 


«  (a)  *  u  0  «■“ 


(4.6) 
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In  this  event 


•lj  r  J“  0 


dej.j  (z') 


(•’) 


The  integral  eq.  (4.7)  is  reducible  to  the  differential  equation 


d*ij  *  270  dz  ®ij  +  270  dsij 


(4.8) 


How,  - — x  ds. .  may  be  identified  as  the  "elastic"  component  of  devjctoric 
2u  0  ij 

strain  of  classical  plasticity.  If  one  follows  the  traditional  definition 

P 


of  "plastic  strain”  de^^  given  below,  i.e., 

deijp  ■  dtij  -  iclj 


(4.9) 


then  in  view  of  eq.  (48) 


deijP  *  2^0  d*  ®ij 


(4.10) 


But  these  are  the  Prandtl-Reuss  relations.  Hence  our  present  theory  contains 
these  relations  as  a  special  cose.  Where  then  does  it  differ  froca  this  theory? 
It  does  in  the  interpretation  of  the  proportionality  coefficient  dz.  In  the 
Prandtl~Reut j  theory  d £  may  be  positive  negative  or  zero,  and  in  fact,  s  has 
bean  identified  with  the  yield  surface,  i.e.,  plastic  action  is  assumed  to 
occur  when  ds  >  0,  where 


s  *  s(s^) 

But  that  dc£j^  is  sero  whenever 


ds  <  s 


(4.11) 


(4.12) 


In  the  present  theory  ds  is  always  positive  if  the  material  is  deforming  (it 
Is  sero  only  when  deformation  does  not  take  place).  Thus  always. 
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dz  >  o 


(4.13) 


Furthermore  dz  is  not  given  by  eq.  (4.11)  i.e.  it  is  not  related  to  some 
yield  surface  but  its  definition  is  entirely  kinematic.  Thus,  no  yield 
phenomenon  or  surface  are  postulated  here.  One  obtains  the  stress  response 
by  merely  monitoring  the  history  of  strain. 

Also  the  theory  admits  a  further  generality  since  u(z)  need  not  consist 
of  a  single  exponential  term. 

For  instance  u  (z)  may  be  of  the  form 

^  u  (z  )  =  y  Q  +  u  1e“etZ  .  (4.14) 

In  this  case,  however,  the  differential  form  of  eq.  (4.7)  becomes: 


2(  u  o  +  v  de^  +  2u  0  a  e^.  dz  =  ds.j  +  a  s^  dz  .  (4.15) 

The  shear  modulus  y,  at  z  =  0,  (initial  modulus),  is  (u  +w  ^  The 
"plastic"  components  of  the  deviatoric  shear  strain  tensor  are  given  from 
eq.  (4.9),  i.e.. 


dz  s 


ij 


(4.16) 


Note  that  eq.  (4.16)  does  not  satisfy  the  Prandtl-Reuc.1?  relations,  which  are 
also  violated  if  one  adds  more  exponential  terms  to  the  right  hand  side  of 
•q.  (4.14),  In  fact  these  relations  will  be  satisfied  if  and  only  if  u  (z  ) 
is  given  by  eq.  (4.6),  i.e.,  u  is  represented  by  a  single  exponential  term 
only.  This  situation  is  not  particularly  disturbing.  Peters  A  Als^**^ 
carried  out  experiments  on  thin  walled  14S-T4  aluminum  alloy  cylinders  by 
loading  these  in  combined  compression  and  torsion  and  found  that  the  Prandtl- 
Reuss  relations  were  not  satisfied,  for  this  particular  metal. 
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Conclusions 

A  theory  has  been  presented  here,  the  scope  of  which  is  wide  enough  to 
allow  a  rational  phenomenological  description  of  mechanical  behavior  of 
materials  under  various  histories  of  strain  and  temperature.  In  particular, 
the  viscoplastic  behavior  of  materials  is  formulated  mathematically,  without 
recourse  to  the  dichotomy  of  the  deformation  history  ir.  plastic  and  clastic 
parts  and  without  the  necessity  of  introducing  discontinuities  in  material 
behavior,  such  as  yield  surfaces. 

The  theory  merely  asserts  that,  to  every  history  of  deformation  gradient 
and  temperature  of  a  neighborhood  there  corresponds  a  unique  state  of  stress 
in  that  neighborhood.  An  entirely  novel  feature  of  the  theory  is  that  these 
histories  are  defined  with  respect  to  a  time  scale,  which  itself  is  a  mate¬ 
rial  property. 

In  this  paper,  we  have  merely  presented  the  framework  of  the  theory 
without  actually  evaluating  the  material  functions  involved,  through  the  use 
of  experimental  data.  This,  however,  will  be  done  in  Part  II  of  this  paper, 
where  it  will  be  shown  that  the  theory  describes  experimentally  observed 
plastic  behavior  of  metals  with  a  remarkable  degree  of  accuracy. 

Stipulation 

This  manuscript  is  submitted  for  publication  with  the  understanding  that 
the  United  States  Government  is  authorized  to  reproduce  and  distribute  re¬ 
prints  for  governmental  purposes. 
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Appendix  I 

The  following  is  a  short  account  of  the  experimental  work  on  (a)  the 
effect  of  the  definition  of  yield  on  the  shape  of  the  yield  surface  and  (b) 
of  the  work  on  viscoplasticity.  The  references  given  are  b>  no  means  exhaus¬ 
tive  and  the  author  wishes  to  apologize  to  people  of  whose  work  he  is  not 
currently  aware. 

In  Ref . ’s  6  and  7,  aluminum  alloy  tubes**  were  subjected  to  shear  pre¬ 
strain  by  twisting  well  into  the  plastic  region  by  a  predetermined  amount. 

The  yield  surface  corresponding  to  this  degree  of  prestrain  was  established 
by  loading  the  tubes  in  combined  tension  and  torsion. 

In  Ref.  6,  Naghdi  fc  ;nd  that  subsequent  yield  surfaces  distorted  in  the 

direction  of  the  shear  axis  with  a  pronounced  Bauschinger  effect  in  shear 

* 

but  there  was  no  effect  on  the  yield  stress  i»i  tension  (i.e.  the  yield  locus 
did  not  change  in  the  vicinity  of  zero  shear  stress). 

On  the  other  hand,  in  Ref.  7,  Ivey  observed  that  the  yield  surface,  in 
addition  to  distortion,  underwent  a  large  amount  of  translation  in  the  direc¬ 
tion  of  the  shear  axis,  so  that  for  large  prestrains  the  origin  of  the  stress 
space  was  outside  the  yield  surface.  However,  he  was  in  agreement  with 
Naghdi  in  that  the  presence  of  shear  prestrain  did  not  affect  the  yield  stress 
in  tension.  Both  authors  used  deviation  from  linearity  in  a  stress  strain 

diagram  as  a  definition  of  yield. 

(8) 

Mair  and  Pugh  to  check  the  absence  of  "cross  effect" ,  carried  out 
their  own  experiments  on  copper  with  a  high  degree  of  isotropy.  However  they 
used  a  different  definition  of  yield,  this  being  the  point  of  intersection 
of  the  initial  straight  part  of  the  stress-strain  curve  with  a  backward  linear 
extrapolation  of  the  "plastic"  part  of  the  stress-strain  curve. 
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Their  results  varied  significantly  from  those  of  Ivey  and  Kaghdi.  They 

found  that  expansion  and  distortion  of  the  initial  locus  took  place  with  a 

strong  cross  effect  between  shear  and  tension.  Also  a  pronounced  Eauschinger 

effect  in  torsion  was  found  with  large  initial  positive  pretorsion.  These 

authors  also  observed  pronounced  "plastic"  unloading  in  shear. 

(9) 

The  results  of  Szezepinski  and  Kiastkowski  tend  to  confirm  the  find- 
(8) 

ings  of  Mair  and  Pugh'  .  Their  results,  moreover,  were  significant  in  other 
respects.  Specifically,  using  the  proof  strain  to  define  yield,  they  studied 
aluminum  alloy  sheets  under  biaxial  tension  with  the  intention  of  finding 
the  effect  of  prestrain  on  the  shape  of  the  yield  surfaces.  They  observed, 
migration,  distortion,  expansion  and  sometimes  rotation  of  the  initial  yield 
locus. 

Similar  conclusions*  can  be  drawn  from  Szezepinski 's  paper as  well 

(17) 

as  Kiastkowski  and  Szezepinski 's  ,  in  which  tubular  brass  specimens  were 
subjected  to  combined  axial  and  circumferential  stress. 

Initial  and  subsequent  yield  loci  were  plotted  when  yield  was  defined 
(a)  as  departure  from  linearity  or  (b)  when  it  was  set  to  correspond  to  a 
certain  proof  strain.  In  particular,  when  definition  (a)  was  used,  subse¬ 
quent  yield  loci  did  not  contain  the  initial  locus,  but  when  (b)  was  used, 
with  ptvof  strain  set  at  0,5%  subsequent  loci  c'--*  "  ’  the  initial  locus. 

Attempts  to  describe  the  change  of  the  yield  locus  with  prestrain,  by 

(13) 

simple  models  have  not  proved  satisfactory.  Batdorf  and  Budianski  , 
suggested  that  after  prestrain,  the  yield  locus  is  the  minimum  surface 

*  In  this  connection,  see  also  work  of  the  samo  general  nature  by 
Bertsch  and  Findley and  Hu  and  Bratt^12\ 
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through  tho  point  of  prestrain  and  the  initial  yield  locus.  This  model 

however  does  not  account  for  the  Bauschinger  effect.  The  kinematic  hardening 

(14) 

rule**  proposed  by  Prager  ,  was  partially  successful,  in  so  far  as  it  can 
be  of  value  only  when  the  stress-strain  curve  of  a  material  in  simple  tension 
is  bilinear^^.  Otherwise  subsequent  shapes  of  yield  locus  must  be  defined 
in  terms  of  a  parameter  that  depends  on  the  history  of  strain^^  to  obtain 
realistic  unloading  behavior. 

A  more  realistic  model  is  the  one  by  Hodge'  '  which  includes  transla¬ 
tion,  expansion  and  distorsion  of  the  yield  surface.  This  model  covers  all 
contingencies  but  does  not  include  the  history  of  stress  cn  the  shape  and 
position  of  the  yield  locus. 

However,  every  definition  of  yield  gives  rise  to  a  different  yield  sur¬ 
face.  If  we  insist  that  the  increment  of  plastic  strain  is  to  be  normal  to 
the  yield  surface,  then,  for  a  complex  but  specific  loading  history,  each 
such  definition  will  give  rise  to  a  different  plastic  strain  history.  Only 
one  of  these  can  be  the  correct  one. 

So  it  appears  that  through  Eisenberg's  and  Phillip* s^^  mathematical 
description  of  a  yield  surface  has  been  most  promising,  we  must  be  prepared 
to  question,  if  necessary,  whether  the  concept  of  yield  point  and  yield 
surface  are  the  only  way  by  which  plastic  effects  may  be  described,  especially 
in  view  of  the  fact  that  these  may  take  place  immediately  following  the  ini¬ 
tiation  of  deformation  of  material,  though  they  may  be  negligible  in  the 
region  of  small  strains.  This  would  agrte  with  the  point  of  view  that  dis¬ 
locations  (and,  therefore,  plastic  behavior)  originate  immediately  upon 
initiation  of  the  loading. 


**  See  also  Ref.  (19) 
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Viscoplasticity 

The  need  for  the  development  of  the  theory  of  viscoplasticity  arises 
from  the  recognition  of  the  strain  rate  sensitivity  of  metals  under  dynamic 
loading. 

The  difficulty  in  trying  to  synthesize  a  rational  "rate"  theory  from 
experimental  observations,  a  priori,  lies  in  the  fact  that  under  dynanic 
conditions  the  inertia  effects  are  significant.  In  the  absence  of  a  consti¬ 
tutive  theory,  these  effects  cannot  be  calculated.*  Therefore,  in  the  case 
of  dynamic  theories,  such  as  viscoplasticity,  theory  and  experiment  must 
advance  together. 

The  literature  abounds  with  data  on  the  subject  of  strain  rate  sensiti- 

(20-30)  (32) 

vity,  particularly  in  one  dimension.  Lindholn  carried  out  dynamic 

% 

experiments  in  one  and  two  dimensions  in  an  attempt  to  generalize  results 

which  were  arrived  at,  by  consideration  of  thermally  activated  processes  and 

their  relation  to  dislocation  theory  in  metals.  See  Ref.'s  35-40. 

An  early  attempt  at  a  theoretical  viscoplastic  constitutive  equation  in 

(31  32) 

one  dimension  is  due  to  Malvern.  *  This  equation  assumes  the  existence 

of  a  "static"  stress-strain  relation  and  then  relates  the  stress  increment, 

with  respect  the  static  value,  to  the  strain  rate. 

Modifications  and  generalizations  of  Malvern's  equation  were  made  by 
(33) 

IiUbliner  who  included  a  limiting  maximum  stress-strain  curve,  and  by 
(34)  (35) 

Persyna  ,  Perzyna  and  Wojno  who  proposed  a  multiaxial  generalization 

*  Constant  strain  rate  experiments  would  appear  to  be  an  exception,  by 
being  less  susceptible  to  inertia  effects.  However,  Ref.  38  tends  to  negate 
this.  Long  specimens  give  different  responses  to  short  ernes,  under  the  same 
conditions* 


for  finite  strains,  assuming  the  additivity  of  the  elastic  and  plastic  strain 

( 39 ) 

components  and  by  Perzyna  who  used  concepts  of  internal  coordinates  and 
irreversible  thermodynamics  to  eliminate  the  above  assumptions  and  to  put  the 
theory  on  sounder  foundations. 

(39) 

Though,  in  his  last  treatment,  Perzyna  abandoned  the  additivity  of 
plastic  and  elastic  strains,  he  still  retained  the  concepts  of  yield  stress 
(and  yield  surface)  and  the  hypothesis  of  a  datum  plastic  stress  strain  rela¬ 
tion,  with  respect  to  which  "strain  history"  is  to  be  related  to  the  "excess" 
stress  through  the  internal  coordinates. 

Qyr  theory  differs  from  Perzyna’ s  theory  in  this  respect. 

We  close  by  mentioning  that,  with  the  exception  of  the  papers  by 
(39) 

Perzyna  ,  only  a  moderate  research  effort  has  been  made  in  the  area  of 
coupling  between  a  viscoplastic  and  a  thermal  process.  However,  Chidister 
and  Malvern^ 25 ^ ,  Lindholn/27^  and  Trozera,  Sherby  and  Dorn ^ 37 ^ ,  considered 
the  effect  of  a  change  in  uniform  temperature  on  viscoplastic  behavior,  with 
9.  view  to  confirming  some  results  of  the  dislocation  theory. 
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sendix 


Thn  form  of  eq.  (3.22)  for  isotropic  materials,  according  to  eq. 's  (3.22 


a-e),  is 


*  C  %  A1  CiiEjj  *  **  A2  £ijEij  +  BlEiiqij  +  B2  Eij  q' 


11  il 


,  _oo  a  .  _a  a  a  .  -a  _a©.  a 

’  18  C1  qii  Ojj  +  55  C2  «ij  ”>j  ♦  d9*M  *  E  Sqii 


+  h  r92 


(A. 2.1) 


As  a  result  eq.'s  (3.8)  and  (3.9)  yield: 


°ij  =  A1  *ij  ‘kk  *  A2Cij  *  B1  Cij  qkk  *  B2  qij 
tDS‘ii  (o  summed) 

*X  *  B  +  E°q®^  ♦  1*9  (o  summed) 


(A.2.2) 


(A.2.3) 


On  the  other  hand. 


*  **  *»  *lj  *  B2  ‘lj  4  C“  *ij  \k 


♦  Cj  q8^  +  E8  (  not  summed) 


(A.2.3) 


Hence,  use  of  eq.  (3.21)  in  accordance  with  eq.  (3.23f)  yields  a  set  of  first 
order  differential  equations  for  q8^ ;  these  can  be  expressed  as  a  set  for 
qJJh  and  another  for  the  deviatoric  part  of  q8^,  which  we  denote  by  p^.  Thus 
in  the  notation  of  eq.'s  (3.31)  and  (3.33) 


*o  ekk  +  Co  qkk  ♦  F*  *  bo  SfiT  *  0 


a®  *  *  <*  .  v a  dplj  , 

82  *ij  +  C2  *ij  4  b2  dT^*  *  0 


(A. 2.4) 


(A.2.5) 
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In  both  eq.'s  (A. 2. 4)  and  (A. 2. 5)  a  is  not  summed.  It  follows  from  the 


abov(e  two  equations  that 


B  rz 


=  -  —  f  5Xo(z‘z,)  r  (z»)dz*  -  ~  f  eXe(z”z,)3(z')dz*  (A. 2. 

.a  I  kk  .  o  I 

b  Jz  b  *z 

~  0  0 


0  0 
_a 

B„  tz 


..-.if  ;p«(z-!:')e..(z')d2- 


(A. 2. 7) 


J  b“  Jz  J 

2  0 

where  X  and  p  are  given  by  eq.  (3.32). 

8  O 

In  the  light  of  the  tensorial  notation  that  we  have  adopted  for  the 


internal  variables,  let 


where,  ||e.,.||  s  |e.^  Cj.|  ,  etc.  Evidently  , 

A2  s  hj  +  1  fi2 

1  3  0 

Then  as  a  result  of  eg.s  (A. 2. 7)  and  (A. 2. 8) 


I  &\W$H 


(A. 2. 8) 


(A. 2.9) 


(A.2.1C) 


' v  .  .  'r 


0  c0 

where  as  beforo  *6 

!  l««P 

Also  from  eq.  (A.2.S), 


*  1  *u  i 


(A. 2*11) 


(A.2.12) 


However,  since 


’ii  *ii  I  ilk  II  Ik  II 


(A.2.13) 
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it  follows  from  (A. 2. 12)  that 


1 


2a 

dr 


i  2Ai 


(A. 2. 14) 


Also  as  a  result  of  eq.  (A. 2.4) 


ft  i 2 1°: 


A  +  2  E  6 
0 


(A. 2. 15) 


At  this  point  we  order  our  internal  variables  as  shown. 


Pij  »  Pij . pij  5  qkk  *  qkk  *  *  *  '  qkk  * 

Let  qa  be  a  typical  internal  variable.  Then,  whether  it  belongs  to  the  p- 

group  or  the  q-group  above,  as  a  result  of  eq.'s  (A. 2. 10),  (A. 2. 11),  (A.2.14) 

and  (A.2.14),  given  two  positive  members  6^  and  6^,  however  small,  we  can 

choose  A  and  A,  (and  therefore  A)  and  A  such  that 
0  A 

ki  i  h and  |ar  |  -  *2  * 
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ABSTRACT 

The  endochronic- theory  of  viscoplasticity  developed  previously 
by  the  author  is  used  to  give  quantitative  analytical  predictions  on 
the  mechanical  response  of  aluminum  and  copper  under  conditions  of 
complex  strain  histories,  One  single  constitutive  equation  describes 
with  remarkable  accuracy  and  ease  of  calculation  diverse  phenomena,  such 
as  cross-hardening,  loading  and  unloading  loops,  cyclic  hardening  as 
well  as  behavior  in  tension  in  the  presence  of  a  shearing  stress,  which 
have  been  observed  experimentally  by  four  different  authors, 


1 ,  Introduction 


In  Part  I  of  this  work,  a  theory  of  viscoplasticity  (of  which  the 

theory  of  plasticity  was  a  part)  was  developed  on  the  basis  of  the  concept 

that  the  current  state  of  stress  is  a  functional  of  the  entire  history 

of  deformation  and  temperature, 

but  the  history  was  defined  with  respect  to  a  time 
.  scale  which  is  in  itself  a  property  of  the  material 
at  hand, 

In  particular  for  the  case  of  strain-rate -independent  materials  under 
isothermal  conditions,  which  is  the  topic  of  this  part  of  our  work,  it 
was  shown  that  within  the  restriction  of  small  strains 

f*  ^"kk  f 

°ij  *  6ij  J  ^z“z’ ^  dz'  +  2  J  y(z-z')  dz'  (1,1) 


where 


1U>  .1  *!  X  .-Or* 

r=l  r 

n  •*ot  * 

y(z)  =  y  +  t  u  e  r 


(1.2) 


(1.3) 


where  Aw,  Ap,  ur»  br  and  are  positive  constants  and 


QZ 

z*z(Oj  ^  >  0,  z  >.0, 

The  symbol  z  denotes  a  positive  monotonically  increasing  time  scale 
with  respect  to  a  time  measure  d;  such  that 


(1.4  a,b) 


*  Pijkt  d‘lj  dtkl  <l-S> 

where  p^*  i*  a  material  tensor,  which  is  positive  definite  and  which,  for 


1 


the  isotropic  materials  envisioned  in  Eq,  (1,1),  has  the  form 

.  Pijkt  S  kl  6ij  6kt  +  k2  6ik  6it  (1,6) 

k2 

where  and  k2  are  material  constants  ,  such  that  k^  +  >_  0,  k2  >_  0. 

It  is  evident  from  Eq,  (1,5)  that  d5  is  independent  of  the  natural 
time  scale  given  by  a  clock  and  thus  materials  described  by  Eq,  (1,1)  are 
strain  history  dependent  but  strain-rate  independent. 

The  derivation  of  constitutive  equation  (1.1)  was  given  in  detail  and 
its  relation  to  classical  theory  of  plasticity  was  examined,  in  some  of  its 
aspects,  in  Part  I« 

the  present  paper  we  shall  be  concerned  with  the  reeil  behavior  of 
metals  under  conditions  of  room  temperature  and  slow  straining.  By  exam¬ 
ining  data  on  copper  and  aluminum  which  were  obtained  in  the  laboratory 
by  various  experimenters,  we  shall  show  that  Eq.  (1,1)  does  indeed  have  the 
capability  of  explaining  quantitatively  and  with  remarkable  accuracy  such 
diverse  phenomena  as  cross -hardening  i.o,  hardening  in  tension  due  to  torsion, 
loading-unloading  loops,  and  hysteresis  loops  during  repetitive  tension¬ 
unloading  -  compression-unloading  histories,  as  well  as  behavior  in  tension 
in  the  presence  of  shear  stress  due  to  torsion, 

2,  Discussion  of  Equation  (1,1) 

At  the  present  time  one  cannot  find  sufficient  experimental  data  in 
the  literature  to  determine  the  functions  z(C),  f(z)  and  u(x).  Therefore  in 
order  to  use  the  theory  at  all,  with  the  object  of  interpreting  available 
experimental  data,  we  have  to  seek  other  avenues,  essentially  heuristic,  to 
determine  the  form  of  the  above  functions. 
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(2.7a) 


then, 


o.2  »  a  oa  a 


'■  <V  ♦  -r  >  «u*5jj 


(a  summed) 


dYy  0  dq., 

eo  3t  *  (bi  *  T,>  f<t>  “ - dib  1 

(a  summed) 

At  this .point  various  possibilities  present  themselves.  The  simplest  is 
to  take  f(c)  =  constant,  in  which  case, 


(2.7b) 


lS{l?+io 


(2.7c) 


where  (  and  c,  are  constants.  However  it  can  be  shown  (and  will  be 
o  1 

demonstrated  in  later  sections)  that  this  choice  eliminates  "cross-hardening" 
in  the  sense  that  a  change  in  the  uniaxial  stress-strain  behavior  due  to 
shear  prestrain  cannot  be  accommodated.  This  effect  has  been  observed  *10 
practically  all  experiments  on  metals  that  have  ever  been  reported.  There¬ 
fore,  though  Eq,  (2,7)  is  convenient,  it  is  not  very  useful, 

This  observation  applies  to  Ilyushin's  and  Rivlin’s  theory,  where 

s  *  ;  ♦  C  and 
o 


r5  i 

C*  J  (de^de^)' 


<2.7d) 


Mote  that  Eq,  (2.7d)is  a  particular  case  of  Eq.  (2.6))  also  note  the 


absence  of  the  material  constants  k^  and  kg  from  Eq,  (2,%)iwhich  renders 
C  independent  of  the  material  at  hand.  In  other  words  their  theories  are 
not  endochronic, 

The  next  natural  choice  ia  to  consider  f({)  to  be  a  linear  function 
of  {  i.e. 


f(t>  «  1  ♦  BC 


(2.6) 


where  8  is  a  positive  constant.  Note  that  8  >  0  because  b  0  >  0,  as  well 

% 

as  -g-r  >  0,  thus  necessitating  that  f(c)  >  0,  for  all  C.  As  a  result  of 
Eq.'s  (2.7)  and  (2.8), 


z  «  j-  log  (1  +  8c) 


(2.9) 


an  expression  which  has  been  found  to  give  excellent  agreement  in  the  cases 
of  soma  significant  experiments,  as  will  be  shown  in  subsequent  Sections. 

In  the  absence  of  experimental  data,  the  question  of  the  form  of  the 

"relaxation"  functions  X(z)  and  u(z)  is  equally  difficult. 

There  are  two  simplifying  assumptions,  however,  which  lead  to  a 

relation  between  X(z)  and  y( z),  so  that  one  is  left  with  the  problem  of 

finding  the  form  of  only  one  of  these  functions.  One  is  that  of  an  elastic 

% 

hydrostatic  response  and  the  other  is  the  assumption  of  constant  Poisson's 
ratio. 

Efficient  use  of  the  first  assumption  is  made  by  writing  Eq.  (1,1) 
in  terms  of  the  hydrostatic  and  deviatonic  components  of  o^,  in  which  case 


Jkk  s  3  f  K(*-*')  *Ckk  dz' 

Jo  V#'  ' 

*  2  J  u(z-t')  dz' 


(2.10) 


(2.11) 


where  K(i)  is  the  bulk  modulus.  Elastic  hydrostatic  response  implies  that 
Hz)  *  Hl(z),  in  which  case  Eq,  (2,10)  becomes, 


"kk  *  3l<Ekk 


(2.12) 


The  assumption  of  constant  Ebisson  ratio  leads  to  the  conclusion  that 
y(x)  and  K(x)  differ  by  a  multiplicative  constant,  and  can  both  be  written 
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in  terms  of  a  single  function  G(z),  such  that 

K(z)  «  K  G(z)  (2.13) 

o 

g(z)  =  U  G(z)  (2.1H) 

o 

1  ,  where  G(0)  =  1, 

This  assumption  has  the  added  advantage  that,  under  condition  of 
plane  stress,  or  uniaxial  strain,  the  strain  in  the  unstressed  directi  n 
is  related  to  the  strains  in  the  stressed  directions  by  a  multiplicative 
constant.  Thus  the  strain  increments  in  the  direction  of  zero  stress 
may  be  easily  eliminated  from  the  express ion  for  dc  so  that  the  latter  may 
be  expressed  solely  in  terms  of  the  strain  increments  in  the  stressed 
directions. 

3 .  CrosshardeninF-_ip_  T_englon=lfflcaiaiL 

It  has  been  observed  that  in  aluminum  and  copper  as  well  as  in  other 
metals,  prestrsining  in  torsion,  well  into  the  plastic  range,  has  a  signi¬ 
ficant  hardening  effect  on  the  stress  strain  curve  in  tension, 

In  this  Section  we  shall  analyze  data  by  Hair  and  Pugh,  who  have 
investigated  this  effect  on  annealed  copper.  Their  experiments  were  per¬ 
formed  accurately  and  with  care,  on  '«ry  thin  circular  cylinders  which 
were  twisted  well  into  the  plastic  region,  so  that  upon  unloading  there 
remained  a  permanent  residual  shear  strain.  The  effect  of  initial  shear 
prestrain  on  the  tensile  response  was  then  obtained  by  loading  the  cylinders 
in  tension. 

The  constitutive  equations  pertinent  to  the  above  situation  are  easily 
found  to  be: 
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0  S 


/; 


E(z-z* ) 


de 

a? 


dz' 


(3.1) 


t  «  2  f  *  M(z-a')  §£,  dz'  (3.2) 

y  o  ^ 

where  o  and  c  are  the  axial  stress  and  strain,  respectively,  and  t  and 
n  are  the  respective  shear  stress  and  tcnsional  shear  strain;  the  moduli 
E(z)  and  p(z)  are  interrelated  through  the  bulk  modulus  K(z).  Their  re¬ 
lating  is  best  expressed  through  their  Laplace  transforms: 


3fc 


(3.3) 


To  deal  with  the  effect  of  cross-hardening  analytically,  we  have 
assumed  a  constant  poisson  ratio.  As  a  result  Eq,  (3.3)  reduces  to  the  form: 


E(z)  *  E  G(z)  (3.4) 

o 

where 


Regarding  the  fora  of  G(s)  we  have  taken  the  simplest  possible  view 
by  assuming  that 

GU)  ■  •*“  (3,6) 

Despite  these  simplification*  we  have  been  able  to  obtain  excellent 
agreement  with  experimental  data  that  have  hitherto  lacked  analytical 
representation. 

Analysis 

In  the  tens ion- tors ion  test  the  effect  of  constant  poisson  ratio  is  to 
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(3.7) 


2 

reduce  d;  to  the  form 

dc2  *  k.  de2  ♦  k2  d  2 
x  2  n 

where  k^  and  k2  are  material  constants,  not  v.,o  same  as  those  in  Eq,  (1.6), 
During  torsion  (c=0), 


C  =  k2n 

wheroas  during  tension  (n=no)  and  after  pretension 


(3,8) 


<  >  k2ne  4  kxc 

where  nQ  is  the  maximum  shear  prestrain. 

Equation  (3.1)  may  now  be  written  in  the  form 


o  *  E 
o 


(3.9) 


(3.10) 


(where  C  *  kn_)  when  allowance  is  made  of  the  fact  that  t  *  0  In  the 

O  2  0 

range  0  <  C  <  k  i  *  Thus  cross -hardening  is  taken  fully  into  account  by 
Eq.  (3,10),  through  the  shear  prestrain  parameter  CQ,  which  appears  as  a 
lower  limit  on  the  integral  on  the  right  hand  side  of  Eq.  (3.10),  If, 
in  particular,  we  assume  that  G(C)  is  given  by  Eq*  (3,6)  and  use  of  this 
is  made  in  Eq,  (3,10)  the  latter  becomes 


o  *  E  a 
o 


-os(c) 


c 


.«<«’)  If. 


(3.11) 


The  integral  in  the  right  hand  side  of  Eq,  (3,11)  can  be  evaluated 
•illicitly  by  using  Eq.  (2,9)  and  noting  that  during  monotonically  increasing 
txwnsion  *  k^  .  Omitting  the  algebra. 


(3.12) 


E  <1+00  .  .  -n 

0  *  THE” 1 11  '  (rfsf  ’  1 

JL  O 


whore 


D  =  1  + 


T 


(3.13) 


and 

C  <  C  «  «  ♦  (3.14) 

Equation  (3.12)  represents  a  family  of  stress  strain  waves  in  tension, 
in  terras  of  the  prestrain  parameter  and  the  “cross-hardcnins”  parameter*  8. 

To  determine  the  material  parameters  in  Eq,  (3.12)  we  note  that  in  the 
absence  of  shear  prestrain  (CQ  *0), 


E  (1+B  c) 

0  s  -  -  - - {  1  -  (1+0  c) 

p  n  1 


(3.15) 


where  6^  *  k^B, 

It  may  be  verified  that  as  e  0,  o  =  E^t  i.e.  Tq  is  the  initial  shape 
of  the  stress-strain  curve.  Also  as  c  increases,  o  tends  syraptotically  to 
the  linear  expression 

0  >  1:0  (its,.)  (3.16) 

1 

*  Thera  is  ample  justification  for  calling  0  the  cross-hardening  parameter, 
Indeed  in  the  limit  of  8  *  0,  and  using  Eq.  (3.9)  Eq,  (3,12)  becomes; 


o  .  pi.  (l-e*8*) 


¥ 


which  is  independent  of  in  other  words  cross -hardening  cannot  taka 
place  when  8  *  0,  as  pointed  out  earlier. 
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It  is  the  slope  (tangent  modulus)  of  the  asymptotic  straight  line, 


then 

E 

o 


(3.17) 


Also,  as  shown  in  Fig.  1,  if  one  extrapolates  backwards  the  asymptotic 
straight  line  to  intersect  the  stress  axis  one  obtains  an  intercept 
from  which  is  determined  by  the  relation 

K 

8,  *=  •—  (3.18) 

1  °o 

Similarly  integration  of  Eq,  (3,2)  yields  an  equation  analogous  to 
Eq.  (3,16) j  this  is 


t  = 


2uo(i+62n 

~ 7~~ 


{ 1  -  (  02nrn) 


(3.19) 


where  &2  *  k20,  Thus,  S2  and  J4q  may  be  determined  from  Eq,  (3,19). 

Finally  we  observe  from  Eq,  (3.12)  that  the  intercepts  '  in  the 
presence  of  shear  prestrain  are  given  from  the  expression 


o  '  =  o  (1+BC  )  =  o  (1+B_n  ) 
o  o  o  o  2  o 


(3.20) 


Equation  (3,20)  was  used  to  confirm  the  self-consistency  of  the  theory, 

*1 

However  Eq.'s  (3,12),  (3,1  )  and  (3,19)  can  only  yield  the  ratio  («-—) 

2 

but  the  constants  k^  and  k^  cannot  be  evaluated.  In  this  sense,  and  for 
these  experiments  one  may  choose  k2  arbitrarily;  we  chose  k2  =  1. 

Experimental  data  obtained  by  Mair  and  Pugh  that  illustrate  the  effect 
of  cross -hardening  are  given  in  Fig.  2, 

Curve  0  is  the  virgin  stress  strain  curve  for  the  type  of  copper  they 
used.  The  circles  on  the  curves  A,  B  and  C  are  experimental  points  corres- 


-2  -2  -2 

ponding  to  initial  shear  prectrains  of  ,25x10  ,  1,5x10  and  3x10 

respectively. 

From  curve  0,  Eq  =  lHxlO6  lb/ In2,  g,  =  ,53xl02,  n  =  M6,  With 
=  l,  £q,  (3,20)  vas  used  to  give  k  =  1,00,  The  curves  A,  B  and  C 
were  then  calculated  and  plotted  as  shown.  Without  a  doubt  the  agree¬ 
ment  between  theory  and  experiment  is  remarkable, 

*t.  Repetitive  loading-unloading  cycles 

The  tensile  strain  history  c(?)  corresponding  to  a  typical  tousile 
loading-unloading  sequence  is  shown  in  Fig.  3,  We  use  the  terms  "straining" 
and  "unstraining"  in  the  following  sense: 

The  ranges  0  <_  c  <  ^  ~  C  *  ;3»  5B  -  C  <  C5  CD  ~  C  rePresent  straining 

in  tension. 

The  ranges  $  <_  ; <  C2.C3  <_  Z  <  tA,  ?5  <_  C  <  represents  unstraining  in 
tension. 

The  ranges  c  <_  C  <  c4*  C  <  Cg,  C^C<  represent  straining  in 
compression. 

The  ranges  54  <.  5  <  ?B,  Cg  <,  C  <  Cy»  <_  ?  <  Cp  represent  unstraining  in 
compression. 

Points  on  the  C-axis  denoted  by  (r-1,2,,,,)  represent  points  of 
discontinuity  in  ^E,  brought  about  by  reverting  from  straining  to  unstraining 

3T 

histories,  or  vice-versa, 

A  perusal  of  experimental  data  on  copper,  shows  that  the  constitutive 
equation  of  the  metal  varies  depending  on  its  previous  history  of  manufacture 

and  subsequent  annealing.  The  single  term  form  of  G(c)  that  explained  Mair 

(2)  (3) 

and  Pugh's  data  so  well  was  found  inadequate  to  explain  data  by  Lubahn 
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(4) 

and  by  Wadsworth  , 

We  found  however,  that  the  adoption  of  a  single  extra  term  for 
G(z)  sufficies  to  describe  quantitavely  broad  trends  of  their  data.  In 
effect  we  took 

G(z)  =  Ej_  +  G2e“az  (4.1) 

or 

E(z)  «  Ex  +  E2e“°Z  (4.2) 

Let  ?n  (m  =  1,2....)  be  the  last  point  of  discontinuity  in  Then  using 

Eq's.  (2.9),  (3.11)  and  (4.1)  and  in  the  range  C  <  ? 

m  "» 

o  =  (1+85)0 „/(-l)“  -  - +  2  J,  (-l)r+1  +  1+E.C  (4.3) 

°\  (1+35)”  r  1  lt6C  J  1 


The  quantities  may  be  evaluated  explicitly  in  terms  of  (the  values 
of  strain  corresponding  to  C  )  by  the  formula 

c  =  2k  J  (-1)8"1  c  +  k.(-l)r  e  (4.4) 

r  1  s-l  si  r 


The  effect  of  E^on  the  unstraining  characteristics  ia  remarkable, 
especially  since  its  effect  on  the  shape  of  the  straining  part  of  the 
stress-strain  curve  is  minimal.  Let  the  history  e(z)  be  one  of  continuous 
straining.  Then  Eq,  (4,3)  becomes! 


E.d+B.e) 

o  =  E.  c  +  ~= — ~ —  (  l-(l+0.c)“n) 
1  nB,  1 


(4.5) 


From  Eq,  (4,5)  we  obtain  the  following  relations  in  the  notation  of 


Section  3. 


(4,6a) 


unknown  material  constants  E^,  E^,  n,  g^,  It  has  been  found  that  a  fourth 
relation  can  be  obtained  by  considering  the  "unloading"  portion  of  the  stress- 
strain  history. 

Fig.  4.  shows  the  stress-strain  relation  for  a  uniaxial  specimen 
which  has  been  strained  in  tension  to  a  strain  value  e^;  whereupon  it  is 
unloaded  and  compressed  until  the  final  strain  is  zero. 

The  strain- intrinsic  time  measure  history  e(c)  corresponding  to  the 
above  stress-strain  history  is  also  shown  in  Fig.  5, 

Equation  (4.3)  in  conjunction  with  the  above  history  yields  the 
relation,  at  e  =  Oj 

2(1+8C  )"-l 

o  *  o(l+80  (  -1+ - =~ -  }  (4,7) 

°  (i+sc)n 

If  the  value  of  is  sufficiently  large  (in  the  case  of  copper  this  value 
was  found  to  be  50x10  \  or  so)  then  is  given  very  nearly  by  the 

expression 

o  *  *  -  0(1+28^.)  =  -  o  (1+28  c  )  (4,8) 

o  o  l  0  1  1 

The  constant  8.  can  now  be  obtained  from  Eq.  (4,8)  and  the  constants  E, , 

Eg  and  n  can  be  found  from  Eq’s,  (4,6a*c), 

We  illustrate  the  points  made  in  the  above  discussion  in  Fig,  6  where  stress 
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strain  curves  for  three  different  materials  are  given  when  these  are 
subjected  to  the  same  strain  history  shown  in  Fig.  5. 

The  constants  for  these  materials  are  given  in  the  following  tables 


E1 

E2 

6 

n 

1 

0 

6.14X106 

,4xl02 

25 

2 

,24x10 6 

5,9xl06 

0 

00 

3 

,12xl06 

6,02xl06 

,2xl02 

50 

What  is  remarkable  is  that  changing  6  results  in  these  materials  having 

indistinguishable  stress-strain  curves  during  straining  but  wildly 

differing  ones  during  unstraining. 

In  Fig.  7  we  illustrate  an  attempt  to  predict  analytically  the 

loading-unloading-loading  response  of  copper  in  simple  tension,  The 

(3) 

solid  line  is  an  experimental  curve  obtained  by  Lubahn  for  a  copper 
specimen  which  had  already  undergone  similar  strain  cycles.  We  have  assumed, 
however,  that  these  have  a  negligible  effect  in  the  response  shown  because 
they  occured  sufficiently  for  in  the  distant  "past". 

The  triangular  points  shown,  wore  obtained  theoretically  from  Eq,  (4,3) 
by  assuming  that  the  specimen  was  continually  extended  (without  unstraining) 
until  the  strain  e  =  51,5x10  was  reached.  The  unstraining-straining 
cycle  was  then  applied, 

Desnite  the  fact  that  E(z)  was  approximated  by  two  terms,  as  in 

Eq,  (4,2)  the  agreement  between  theory  and  experiment  is  remarkable.  The 

3  6  3 

constants  employed  were,  oqs  6x10  ,  s  .125x20  ,  6  »  ,02x10  ,  n  e  160. 

In  fact  we  are  not  aware  of  another  instance  where  an  attempt  was  made 
to  describe  such  experimental  data  analytically  by  means  of  one  single 
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constitutive  equation,  In  addition  ve  can  say  with  assurance  that  the 

observed  difference  between  theory  and  observation  can  be  reduced  further 

by  including  more  exponential  terms  in  the  series  representation  for  E(z), 

We  fconclude  this  Section  by  considering  the  effect  of  work  hardening  under 

•  (4 ) 

cyclic  straining.  In  particular  we  shall  examine  the  work  of  Wadsworths 
and  show  that  our  theory  again  provides  an  excellent  analytical  basis  for 
his  results. 

In  this  work  single  copper  crystals  were  tested  under  conditions  of 

uniaxial  cyclic  strain;  The  data  was  presented  in  terms  of  the  resolved 

shear  stress  and  strain  in  the  plane  of  slip. 

Fig.  8  gives  the  first  few  cycles  of  his  straining  program,  in  which 

-3 

a  crystal  was  cycled  under  fixed  limits  of  resolved  shear  strain  of  7x10  . 

The  "peak  stresses"  corresponding  to  the  extreme  values  of  tensile  and 
compressive  strain  increased  monotonically  with  the  number  of  cycles. 

In  Fig.  9  the  values  of  peak  tensile  and  compressive  stresses  have  been 
plotted  by  Wadsworth  against  fdri|,  Ir  is  rather  interesting  that  he 
felt  that  such  a  plot  was  meaningful,  without  further  elaboration  on  this 
point.  Of  course  |dn|,  but  for  a  scalar  factor,  is  our  intrinsic  time 
measure. 

The  history  of  the  resolved  shear  strain  versus  4  is  shown  in  Fig,  10, 

From  this  Figure  it  follows  that  C  =  (2m-l)Ak  ,  Equations  (4,6c)  and  (4,8) 

mi  • 

were  now  utilized  to  find  Bk^,  which  we  denote  by  8^,  and  E^.  It  was  found 

that  8^  s  12,3  and  E^  =  2xl09  At  this  point  n  could  not  be  determined 

cm 

because  the  initial  slope  of  the  stress-strain  curve  corresponding  to  & 
could  not  be  evaluated  accurately. 

However  letting  s  t | t  it  was  found  that  as  m  ■*  Eq,  (4,3) 
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yields  the  asymptotic  expression: 


t  *  njh  A  +  E,  A 
o  1 


(4.9) 


Hence,  from  the  tensile  experimental  curvo  of  Fig.  4  of  Ref.  4  ,  n  could 

be  determined  explicitly  and  was  found  to  be  equal  to  225,  For  this  value 

of  n  the  term  n  was  found  to  be  negligible  for  m  ^  1. 

m 

Thus  for  the  history  in  Fig.  10,  Eq,  (4,3)  gives 


Vi  *  ti+a^awOi]  (-U  .2  r 


m  l+(2r-l)0.A  n  rtl 

t*1)  ♦  v 


(4.10) 


The  above  equation  can  be  simplified  further  for  large  values  of  m.  In 
particular  for  m  >  50,  it  was  found  that  the  series  in  the  bracket  on  the 
right  hand  side  of  Eq.  (4,10)  degenerates  into  the  geometric  series 


r  (  „\r  « 

nil  <’u>  *  ~  1ST 


(4.11) 


where 


fl+(2m-l)81A  \n 
U  =  (l+(2m+ri81A  ] 


(4.12) 


Equation  (4,10)  may  now  be  written  in  an  asymptotic  form  in  terms  of  the 
absolute  value  of  the  shear  stress  as  follows: 


| T  |  *  (  1+B(2m+1)A  k  1-1+  rfvW.A 
m  0  l+u  1 


(4.13) 


For  very  large  values  of  m  (m  »  n)  Eq*  (4,13)  simplifies  further  and 


becomes 


1  1  "eiV 

I’J'-T-  ♦  V 

J- » 


(4.14) 
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Thus, 


limit  I  =  n3,t  A  +  E  A  (4.15) 

I  m  I  lo  1 

m  « 

In  Fig,  11  a  plot  has  been  made  of  the  theoretical  relation  between 
It^l  and  mA  obtained  from  Eq.  (4,13),  The  experimental  points  obtained 
by  Wadsworth  are  also  shown.  The  following  comments  are  in  order.  Though 
our  theory  does  give  values  for  which  are  different  in  tension  from  those 
in  compression,  the  difference  is  not  as  great  as  the  experimental  data 
indie  ite,  and  is  too  small  to  be  plotted  on  the  scale  shown,  However,  the 
theoretical  curve  lies  very  close  to,  and  is  in  fact  bounded  by  the  exper¬ 
imental  points,  which  indicate  a  deviation  between  the  values  of  compressive 
stress  and  those  of  tensile  stress  which  increases  with  m  but  is  never  greater 
than  5.5%, 

This  is  the  first  time  that  a  theory  of  plasticity  has  provided  a  national 
explanation  for  the  phenomena  of  cyclic  hardening, 

5.  Tensile  response  in  the  presence  of  initial  shear  stress 

In  Section  3  we  obtained  a  theoretical  prediction  of  the  effect  of 
prestraj-u  in  torsion  on  the  stress  strain  curve  in  tension.  In  this  Section 
we  shall  examine  theoretically,  in  the  light  of  our  endochronic  theory,  the 
effect  of  initial  constant  prostress  in  torsion  on  the  stress-strain  curve 
in  tension,  To  do  this,  we  have  assumed,  just  as  we  did  in  Section  3,  that 
F.( z)  and  u(z)  are  proportional  to  some  relaxation  function  G(z),  and  further¬ 
more  that  G(z)  consists  of  a  single  exponential  term  i.e,  it  is  given  by 
Eq.  (3.6).  Thus 

E(z)  *  Eq  e"®*  (5.1) 
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In  the  light  of  Eq,  (5,1)  and  bearing  in  mind  Eq,  (2,9),  the  constitutive 
Eq.s  (3,1)  and  (3,2)  can  be  reduced  to  the  differential  equations 


..  dc  _  ao  '  do 
Lo  HT  ■  I+HT  +  Hr 

0  dn  .  ax  dT 

2vo  nr  "  i+sr  +  nr 

where,  as  in  Section  3, 

2  2  2 
dr  s  k  de  +  k2dn 


(5.2) 

(5.3) 


As  mentioned  above  the  test  to  bo  discussed  consists  of  applying  an  initial 

stress  t°  corresponding  to  an  initial  strain  no»  then  keeping  t°  constant, 

a  axial  strain  e  is  applied  and  the  axial  stress  o  is  measured.  The  object 

at  hand  is  to  deduce  from  Eq.'s  (5, 1-5, 4)  the  relation  between  o  and  c, 

(5) 

and  compare  with  tho  experimental  data  obtained  by  Ivey  * 

To  accomplish  this  we  proceed  as  follows.  From  Eq,  (5,4)  it  is  clear 
that  the  axial  straining  process  begins  at  C  *  CQ  where 


C  =  k  n 
o  2  o 


During  this  process  *  0,  bo  that  from  Eq.  (5,3) 

j  _  at0  dc 

dn  8  TT  IT h 


(5.5) 


(5.6) 


Equations  (5,4)  and  (5,6)  now  combine  to  show  that  during  the  axial  st*ainini 
process , 


dc2  *  k^  dc2  ♦  k22  (5-J-) 


w  (!♦$?) 


At  this  point  we  introduce  the  variable  6  such  that 


(5.7) 
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t  *  kx0 


(5.8) 


Also  let  )=  k,  k^B  =  8  and  c  =  (ketT°/2uo).  Then,  in  terms  of  0  and 
as  a  result  of  £q.  (5.7) 


d6 


(i 

l  (1+6^)  ) 


=  de 


(5.8) 


Equation  (5,8)  may  be  integrated  subject  to  the  initial  condition  that-  at 
5  s  k2nQ,  e  *  Oj  or,  8  =  6q  =  kno,  c  *  0. 

Equation  (5,2)  may  now  be  integrated  with  respect  to  8  to  yield 


o  « 


(l+0,8)(<*/8.)  *8 
1  1  o 


J  (l+f^e*)  (o/6r1)yu+8^ 


e')2-c2def  (5.9) 


We  introduce  now  a  change  of  variablo  by  the  relation 


1  +  8^6  s  c  cosh$ 


(5.10) 


whereupon  Eq,  (5,9)  becomes  I 


E  c 

a  «  -P  /  (coshV  -  coshn’V)d8*  (5.11) 

6,  cosh®*’1* 

whore  before,  n  *  1  ♦  jj-  , 

Now, 

f coshnn  dx  «  J cosh  n"2x  dx  ♦  coshn,,ix  ainx  (5.12) 

Since  for  asymptotically  larbe  n  (say  n  >  30),  -  1,  it  follows 

from  Eq,  (5*12)  that  in  this  instance 

(coshnx  -  coahn~2n)  dx  *  i  cosh*1"1*  sinhx  (5..  ) 

U 
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The  result  of  Eq,  (5,13)  can  be  utilized  to  obtain  a  closed  form  solution 
for  o  which  now  becomes, 

E  c  f  cosh*  n-1  ) 

« *  ( *l"h*  ‘  ‘sskt-  >  5lnh,0  ]  •  <5-“> 

Equation  (5.8)  may  also  be  integrated  with  respect  to  *  to  yield 

C  «  §-{F(*)  -  F(*  )  }  (5.15) 

®1  ° 

where 

F(*)  *  sinh*  -  tan”*  (sinh*)  (5.16) 

Thus  a  and  c  are  related  parametrically  through  *  and  +o  such  that 

.  (1+8. kn  ) 

*  *  cosh**1  { - ■■■---  )  (5.17) 

o  c  • 

The  relation  between  o  and  c  has  been  calculated  with  the  following  values 

of  the  constants: 

h  3  2 

k  *  1,  o  *  17,4x10  lb/ln  ,  8  *  20,  n  «  33. 

o  Pj^n 

The  result  was  compared  with  one  of  Ivey's  experiments  in  which  t°  »  14x10* 
lb/m  ,  nQ  *  2.35x10  ,  The  predicted  and  experimental  stress-strain 

responses  compare  very  favorably.  See  Tig.  12, 

Conclusion 

On  the  evidence  of  the  results  presented  above  it  appears  that  the 
endochronic  theory  of  plasticity  can  predict  ' curat ely  the  mechanical 
response  of  metals  under  complex  straining  histories.  The  full  implications 
of  the  theory  will  be  investigated  further  in  our  future  lock. 
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Fig.  10.  History  of  resolved  shear  strain 
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Abstract 


The  hypothesis  of  plastic  incompressibility,  which  is  pivotal  in  the 
development  of  the  classical  theory  of  plasticity  has  been  tested  by  means 
of  a  series  of  simple  tension  tests  on  aluminum,  copper  and  low  carbon  steel. 
The  experimental  measurements  show  conclusively  that  these  metals  are  plas¬ 
tically  compressible,  the  compressibility  increasing  with  straining  in  the 
plastic  region. 

A  new  phenomenon,  termed  lateral  instability,  has  also  been  observed, 
consisting  in  an  abrupt  plastic  flow  in  the  transverse  direction  either  im¬ 
mediately  preceding  or  immediately  following  yield  in  the  longitudinal  direc¬ 
tion,  during  a  simple  tension  test  on  low  carbon  steel.  To  a  less  pronounced 
extent  such  a  phenomenon  has  also  been  observed  in  the  case  of  commercially 
pure  aluminum. 


I.  Introduction 


One  of  the  assumptions  of  the  classical  theory  of  plasticity  is  that 
the  material  remains  plastically  incompressible  while  it  undergoes  plastic 
deformation.  This  implies  that  the  hydrostatic  response  of  the  material 
remains  elastic,  (linearly,  for  small  strains)  even  after  yielding  has 
occurred.  This  assumption  also  leads  to  the  conclusion  that  Poisson's  ratio* 
tends  to  —  as  the  plastic  strain  increases^ 

In  the  present  investigation  the  above  assumption  is  eval,  '’ted  criti¬ 
cally  by  means  of  a  series  of  simple  tension  tests  on  aluminum,  copper  and 
low  carbon  steel.  We  shall  show  in  the  following,  that  cur  experimental 
data  do  not  support  the  hypothesis  of  plastic  incompressibility.  Moreover, 
Poisson's  ratio  does  not  tend  to  ~  as  the  plastic  strain  increases  hut,  on 
the  contrary,  its  value  decreases  and  becomes  smaller  than  its  "elastic" 
value . 

A  newly  observed  form  of  material  instability.  The  above  statements 
on  the  observed  behavior  of  Poisson's  ratio,  are  true  for  copper  and  broadly 
true  for  aluminum  and  mild  steel.  However,  in  the  case  of  mild  steel  (and 
to  a  lesser  extent,  aluminum)  an  unusual  phenomenon,  which  we  shall  call 
LATERAL  INSTABILITY,  was  observed  for  what  is  believed  to  be  the  first  time. 
At  least,  to  our  knowledge  no  such  phenomenon  has  ever  been  reported  in  the 
literature  in  the  past. 


*It  must  be  pointed  out  that,  in  this  reoort,  we  do  not  regard  Poisson's 
ratio  as  some  fundamental  physical  property  or  function,  but  simply  as  the 
ratio  of  the  algebraic  values  of  the  tr^-^verse  and  longitudinal  strains, 
in  a  uniaxial  test.  On  the  other  hand  i,-. cause  the  experiments  reported  here 
were  carried  out  in  the  small  strain  range,  Poisson's  ratio  is  a  convenient 
quantity  to  work  with,  since  it  can  be  used  to  relate  directly  the  uniaxial 
to  the  volumetric  strain;  the  latter  is  an  important  parameter  in  our  inves¬ 
tigation. 
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Lateral  instability  is  defined  by  one  of  the  following  two  observed 
sequences  of  events  during  a  simple  tension  test. 

(a)  As  the  yield  point  of  a  low  carbon  steel  specimen  is  reached, 

the  material  begins  to  flow  longitudinally  (i.e.,  in  the  direction 
of  applied  stress)  and  one  observes  the  usual  increase  in  longitu- 
di.il  strain  while  the  longitudinal  stress  remains  constant. 

However  an  abrupt  occurence  of  plastic  flow  in  the  transverse 
direction  has  also  been  observed,  just  after  longitudinal  flow 
has  ceased  and  strain-hardening  has  begun. 

^  (b)  Abrupt  transverse  f low  occurred  just  before  the  onset  of  longitu¬ 
dinal  flow. 

It  is  significant  that  no  transverse  flow  has  been  observed  during  the 
process  of  longitudinal  flow.  A  more  detailed  discussion  of  this  phenomenon 
will  be  given  in  text. 

In  the  case  of  aluminum  the  onset  of  latoral  plastic  flow  was  not  abrupt, 
but  what  was  observed  amounted  to  a  significant  increase  in  the  transverse 
strain  during  a  relatively  small  increment  in  longitudinal  strain,  A  quanti¬ 
tative  plot  of  these  events  is  shown  in  later  figures. 

Loading- unloading  loops  hove  also  been  observed  in  some  of  the  tests. 

It  will  be  shown  in  later  sections  that  loops  are  also  obtained  in  plots  of 
transverse  strain  versus  longitudinal  strain,  and  hydrostatic  stress  versus 
volumetric  strain. 


3 


1 1 .  Materials  arid  Specimens 

Altogether  three  metallurgically  important  metals  were  tested.  They 

were  commercially  pure  aluminum,  electrolytic  tough  pitch  copper  and  low- 

carbon  steel.  The  aluminum  specimens  were  sheared  from  an  1100-0  aluminum 

sheet,  all  the  copper  specimens  were  cut  out  of  an  electrolytic  tough  pitch 

copper  110  (99,9+%  copper)  bus  bar,  and  the  steel  specimens  were  cut  from  a 

cold  drawn  C  1018  flat.  The  nominal  dimension  of  all  the  aluminum  and  copper 
1  3 

specimens  was  — "  x  — "  x  15"  and  that  of  all  the  low-carbon  steel  specimens 

was  x  ~"x  12".  All  the  specimens  of  the  same  metai  were  considered  to 
8  2 

be  identical  prior  to  the  heat  treatment.  These  specimens  have  been  divided 
into  several  groups,  according  to  the  annealing  temperatures  and  purposes  of 
study.  They  are  listed  in  Table  I. 
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III.  The  Equipment  and  Strain  Measurements 


All  tests  were  conducted  at  room  temperature  on  a  Tinius  Olsen  Lo  Cap 
universal  testing  machine  with  strain  gage  load  cell  weighing  and  equipped 
with  a  Tinius  Olsen  Model  51  electronic  recorder.  The  load  readings  on  the 
dial  indicator  were  calibrated  against  a  Tinius  Olsen  proving  ring  with  the 
electrical  vibrating  reed. 

A  Tinius  Olsen  S-400-2AB  extensometer  was  used  to  measure  the  longitu¬ 
dinal  strains  of  all  the  specimens  except  those  of  group  III  of  aluminum. 

This  is  an  averaging  type  extensometer  with  a  gage  length  of  two  inches. 
The"accuracy  of  the  extensometer  is  ±  0.0001  in/in  as  guaranteed  by  the 

manufacturer  and  was  calibrated  against  the  SR-4  electric  strain  gages  of 

-5 

type  A-7  which  have  the  accuracy  of  2  x  10  in/in.  The  load  vs.  longitu¬ 
dinal  strain  curves  were  plotted  by  the  Model  51  recorder  during  experi¬ 
ments  . 

To  measure  the  transverse  strains  e^,  a  SR-4  electric  strain  gage  of 
type  A-7  was  transversely  put  on  one  side  of  each  specimen.  The  transverse 
strain  was  then  read  dii'ectly  from  a  SR-4  portable  strain  indicator  at  the 
accuracy  of  2  x  10~5  in/in. 

The  strain  measurements  of  the  specimens  of  group  III  of  aluminum  were 
specially  arranged  for  greater  accuracy.  A-7  typo  strain  gages  were  longitu¬ 
dinally  and  transversoly  glued  on  both  sides  of  eaeh  specimen  belonging  to 
this  group,  in  order  to  cancel  out  the  bending  effect.  The  strain  reading 
of  each  gage  was  then  obtained  from  a  SR-4  portable  strain  indicator  with 
the  help  of  a  BLH  Model  225  switching  and  balancing  unit,  finally,  the 
average  value  of  the  readings  of  the  two  longitudinally  arranged  strain 
gages  gives  us  the  longitudinal  strain  cx,  and  the  average  value  of  the 
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readings  of  the  two  transversely  arranged  strain  gages  gives  us  the  transverse 
strain  e  . 

y 


IV,  Description  of  Tests  and  Results 
Experiments  on  Aluminum 

The  first  group  of  aluminum  specimens  was  tested  under  the  "as-received" 

condition.  During  the  tests,  the  Tinius  Olsen  machine  was  operated  at  a 

-4 

constant  rate  of  ci^osshead  separation.  The  rate  was  8  x  10  in/sec,  which 

-5 

corresponds  to  a  strain-rate  of  approximately  5  x  10  in/in/sec  for  the 
specimens  used.  The  recorder  traced  out  the  load  versus  curve  as  readings 
of  load  and  of  transverse  strain  were  taken  simultaneously  at  predetermined 
readings  of  load,  at  small  longitudinal  strains,  and  at  predetermined  read¬ 
ings  of  longitudinal  strains,  when  the  strain  was  large. 

The  stress-strain  curves  are  shown  in  Fig.  1,  in  which  the  curve  for 
specimen  #1  is  identical  to  that  for  specimen  #2.  Fig.  2  shows  the  relations 
between  the  transverse  and  the  longitudinal  strain,  and  Fig.  3  shows  the 
hydrostatic  stress  (x  3)  versus  the  volumetric  strain  curve  for  aluminum  of 
group  I.  A  loading-unloading  loop  was  generated  at  the  end  of  each  test,  and 
we  obtain  loops,  which  is  similar  to  the  hysteresis  loops,  in  the  curves  of 
Figures  2  and  3  (not  shown).  More  observations  and  discussions  of  the  loops 
will  be  presented  later  in  this  report.  The  large  increase  of  -e  at  a  lon¬ 
gitudinal  strain  of  approximately  2.5%  (see  Fig.  2)  corresponds  to  a  case  of 
"mild"  lateral  instability  as  compared  to  the  "strong"  lateral  instability 
observed  in  low-carbon  steel.  We  shall  mabj  more  observations  later  about 
the  lateral  instability  of  aluminum. 

The  second  group  of  aluminum  specimens  was  tested  exactly  in  the  same 
way  as  those  of  group  I,  except  that  the  specimens  were  first  annoaled  at  a 
temperature  of  600°F  for  two  hours  and  then  oven  cooled  to  the  room  tempera¬ 
ture  before  testing.  This  annealing  process  removes  most  of  the  effects  due 
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to  cold  work,  and  the  material  can  thus  be  considered  as  isotropic.  The 
stress-strain  curves  for  specimens  hi  #3  and  hi  ft 4  are  shown  in  Fig.  4,  in 
which  the  two  curves  are  indistinguishable.  We  notice,  however,  that  the 
material  is  considerably  softened  due  to  the  heat  treatment. 

To  see  the  sudden  increase  in  magnitude  of  the  transverse  strain  at  a 
longitudinal  strain  of  approximately  2.2%,  which  has  been  an  observed  fact 
in  the  test  of  group  I ,  this  series  of  tests  was  carried  out  to  the  extent 
that  the  longitudinal  strain  at  the  time  of  termination  of  the  test  was  ap¬ 
proximately  9%.  Here,  again,  we  have  observed  a  sudden  increase  in  -e  at  a 
longitudinal  strain  of  approximately  2.2%.  This  phenomenon  has  been  observed 
both  for  specimens  hi  #3  and  hi  #4. 

The  -e  vs.  e  curves  are  plotted  in  Fig.  5,  and  the  o,,  vs.  e.  ,  curves 
y  x  kk  kk 

are  given  in  Fig.  6.  The  sudden  increase  of  the  transverse  strain  in  Fig.  5 
gives  rise  again  to  the  sudden  decrease  of  the  volumetric  strain  in  Fig.  6 
which  corresponds  to  the  lateral  instability  mentioned  earlier.  In  this 
series  of  test,  the  loopings  of  the  curves  arc  once  more  present  in  all  the 
figures  concerned  (not  shown).  These  loopings  are  due  to  the  loading-unload¬ 
ing  test  conducted  at  the  end'of  each  study. 

It  is  important  to  point  out,  that  aside  of  the  sudden  change  of 
Poisson's  ratio  immediately  following  yielding,  Poisson's  ratio  remains  ap¬ 
proximately  constant  throughout  the  range  tested.  Poisson's  ratio  in  the 
plastic  range  is  defined  in  the  present  investigation  by  the  expression 

Cy 

y  S  «  wdm 

which,  in  fact,  is  also  used  in  elasticity  (see  footnote  in  the  introduction). 

Group  III  of  the  aluminum  specimens  was  tested  quasi-statlcally,  i.e., 
readings  were  taken  after  the  machine  had  been  stopped  and  time  of  5  minutes 
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had  already  elapsed.  By  doing  this,  it  is  believed  that  the  results  obtained 
are  time  independent.  As  it  has  already  been  described  earlier  in  this  re¬ 
port  that  four  strain  gages  were  applied  to  each  specimen  of  this  series  to 
eliminate  the  bending  effect  due  to  the  possible  eccentricity  of  the  loading 
arrangement.  The  comparison  of  the  present  results  (shown  in  Figures  7-12) 
with  those  obtained  earlier  in  groups  I  and  II  of  the  aluminum  specimens  shows 
that  reliable  data  can  be  obtained  without  the  complicated  strain  gage  set¬ 
ups  of  group  III.  In  fact,  the  grips  of  the  Tinius-Olsen  machine  are  of  the 
self-aligning  type,  therefore,  eccentricity  is  to  be  expected  at  its  minimum. 

Loops  are  again  present  in  all  the  figures  concerned  (Figures  7-12).  It 
is  seen  from  Figures  9  and  10  that  Poisson's  ratio  decreases  again  drastically 
immediately  following  yielding. 

The  fact  that  all  the  curves  in  Figures  3,  6,  11  and  12  bend  over  toward 

the  e.  -axis  leads  to  the  conclusion  that  aluminum  is  plastically  compressible, 

-  L  "  ’  *  1 

at  least  when  the  hydrostatic  stress  is  of  a  tensile  nature.  The  assumption  of 
plastic  incompressibility  in  the  classical  theory  of  plasticity  would  load  to  a 
linear  relation  betweon  hydrostatic  stress  and  volumetric  strain;  this  apparently 

is  not  the  case. 

Experiments  on  Copper 

The  experiments  on  group  I  of  the  copper  specimens  were  designed  to  fur¬ 
ther  study  the  loops  which  we  had  observed  in  the  tests  of  aluminum  specimens. 
These  specimens  of  copper  were  tested  under  the  "as-received"  condition.  The 
rate  of  crosshcad  motion  was  kept  constant  (8  x  10  in/sec)  during  this 
series  of  test.  The  loading-unloading  curves  for  specimens  Cu  #1  and  Cu  12 

are  shown  in  Figures  13  and  m,  ard  the  corresponding  -c  versus  c  curves 

y  * 

are  shown  in  Figures  IS  and  16.  We  see  here  again  that  loops  are  definitely 
in  existence  and  that,  after  a  loop  has  been  completed,  the  curve  of  -c^ 
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versus  e  is  virtually  a  continuation  of  the  original  -e  versus  c  curve, 
x  y  x 

The  plots  of  vs.  c  ^  for  these  specimens  are  given  in  Figures  17  and  18. 

The  copper  specimens  of  group  II  were  tested  in  exactly  the  same  way 
as  those  of  the  aluminum  specimens  of  group  II,  except  that  the  annealing 
temperature  was  750°F.  The  results  are  shown  in  Figures  19,  20  and  21.  We 
remark  that  in  contrast  with  the  aluminum  specimens  we  do  not  observe  the 
sudden  increase  in  the  magnitude  of  the  transverse  strain  for  all  the  copper 
specimens  concerned.  As  a  result,  Poisson's  ratio  lr  approximately  a  con¬ 
stant  throughout  the  range  tested. 

Experiments  were  conducted  on  group  III  of  the  copper  specimens  under 
quasi-static  conditions,  i.e.  readings  were  taken  after  the  machine  had  been 
stopped  and  5  minutes  had  elapsed.  The  results  of  this  series  are  given  in 
Figures  22,  23  and  24.  We  again  observe  no  sudden  increase  in  the  magnitude 
of  the  transverse  strain  and  that  Poisson’s  ratio  is  virtually  constant 
throughout  the  range  tested. 

The  tests  on  copper  specimens  of  groups  I,  II  and  III  show  conclusively 
that  copper  is  also  plastically  compressible,  as  it  is  seen  easily  from 
Figures  17,  18,  21  and  24;  otherwise  these  plots  would  have  been  straight 
lines. 

Experiments  on  Low  Carbon  Steel 

The  tame  experimental  procedure  as  described  above  was  followed  when  the 
five  low- carbon  steel  specimens  were  tested.  All  the  five  specimens  were 
annealed,  prior  to  tests,  at  a  temperature  of  1250° F  for  one  hour  and  then 
oven  cooled  to  room  temperature.  Thus,  all  the  specimens  can  be  considered 
to  be  identical  and  the  material  can  be  considered  to  be  isotropic. 

The  tests  were  carried  out  under  quasi-static  conditions,  and  a  typical 
stress-strain  curve  for  these  specimens  is  given  in  Fig.  25.  In  Fig.  26,  the 
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relation  between  the  transverse  strain  and  the  longitudinal  strain  is  plotted 
for  all  the  five  specimens.  It  is  seen  that  aside  of  the  portion  reflecting 
"lavoral  instability"  the  curve  is  almost  a  straight  line.  Of  all  the  speci¬ 
mens  tested  in  this  series,  three  of  them  (#1,  (?2  and  #4)  showed  lateral 
instability  before  the  occurrence  of  the  longitudinal  plastic  flow,  and  follow 
the  path  OABDE  in  Fig.  26.  The  rest  of  the  specimens  (fi 3  and  if 5)  experienced 
lateral  instability  immediately  after  the  longitudinal  plastic  flow  had  been 
accomplished,  i.e.,  lateral  instability  had  occurred  before  the  material 
strain-hardeneu.  In  this  case,  curve  OACDE  of  Fig.  26  is  followed. 

We  conjecture  that  had  the  material  been  homogeneous  as  well  as  isotropic 
then  longitudinal  and  transverse  flow  would  have  taken  place  simultaneously. 
Since  in  practice  this  is  rarely  true,  the  order  in  which  these  flows  occur 
must  then  be  decided  by  slight  differences  in  the  directional  properties  of 
the  specimen. 

To  ensure  the  correctness  of  our  observation  about  the  lateral  instabil¬ 
ity  as  a  material  property  of  low-carbon  steel.  The  transverse  dimension  of 
all  specimens  was  measured  after  the  experiments  by  meanr  of  a  micrometer. 

The  permanent  transverse  strains  were  then  calculated  referring  to  the  orig¬ 
inal  dimensions  of  the  specimens.  The  results  checked  favorably  (with  less 
accuracy,  of  course)  with  those  obtained  by  means  of  electric  strain  gages. 

Two  of  the  five  experiments  (12  and  #4)  were  terminated  right  after  the  occur¬ 
rence  of  the  transverse  plastic  flow  and  before  much  longitudinal  flow  had 
occurred,  in  order  for  us  to  study  the  lateral  instability  more  closely. 
Measurements  by  means  of  both  electric  strain  gages  and  micrometer  showed 
that  large  amount  of  plastic  flow  did  occur  transversely  at  yielding.  If  it 
was  not  for  th*  transverse  plastic  flow,  the  transverse  permanent  strains 
would  have  been  very  small  for  these  two  specimens. 
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The  relation  between  the  hydrostatic  stress  and  the  volumetric  strain 
is  shown  in  Figures  27-31.  The  alphabetical  order  of  letters  in  the  figures 
denotes  successive  states  of  the  material  during  the  tests.  The  last  letter 
in  the  figure  gives  us  the  position  where  the  corresponding  test  was  termi¬ 
nated.  It  is  seen  that  the  lateral  instability  corresponds  to  a  sudden  de¬ 
crease  of  volumetric  strain.  In  the  case  when  lateral  instability  occurs 
before  the  longitudinal  plastic  flow,  the  volumetric  strain  could  become  nega¬ 
tive!  A  very  strange  result! 

However,  a  similar  phenomenon  was  observed  by  Eridgman^^.  He  investi¬ 
gated.,  by  means  of  a  dilatometer,  volume  changes  during  simple  compression 
in  the  plastic  range,  in  various  materials,  such  as  mild  steel,  Norway  iron, 
cast  iron  and  rock  (soapstone,  marble  and  diabase).  He  observed  an  increase 
in  volume  under  conditions  of  compressive  hydi'ostatic  stress,  when  the  axial 
compressive  stress  exceeded  its  yield  value. 
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V.  On  the  Accuracy  of  the  Experimental  Measurements 

It  is  to  be  noted  that  is  is  impossible  to  give  a  strain  gage  on  the 

specimen  such  that  the  strain  gage  makes  an  angle  of  exactly  90°  degrees 

with  the  longitudinal  direction  of  the  specimen.  Slight  errors  are  always 

possible,  and  these  errors  tend  to  make  the  readings  of  transverse  strain 

smaller  than  it  should  be.  These  are,  however,  small  since  the  gages  are 

[3] 

in  the  principal  strain  directions.  According  to  Perry  and  Lissner  ,  a  2 
degree  error  in  gage  alignment  would  only  result  in  an  error  of  less  than 
1%. 

Our  purpose  in  the  present  experimental  study  is  not  the  determination 
of  the  exact  value  of  Poisson's  ratio.  Our  primary  concern  is  the  trend  in 
the  variation  of  Poisson's  r  tio  s  the  longitudinal  strain  increases.  For 
this  reason,  the  above  mentioned  deviations  due  to  the  misalignment  of  the 
transverse  strain  gages  are  not  important,  since  they  hardly  affect  the  broad 
trends  in  the  variation  of  Poisson's  ratio,  which  are  of  interest  here. 

Slight  scattering  of  some  of  our  experimental  data  for  different  speci¬ 
mens  may  be  attributed  partially  to  the  above  mentioned  misalignments  of  the 
transverse  strain  gages.  We  nctice,  however,  that  our  data  are  very  consis¬ 
tent  in  general,  and  the  slight  scattering  of  data  occurs  only  in  the  case 
of  copper  specimens  H 5  and  #6.  Even  in  these  cases  where  data  are  scattered 
slightly,  the  curve  for  each  specimen  is  itself  smooth,  which  indicates  that 
scattering  is  due  to  a  largo  extent  to  variability  in  the  properties  of  the 
specimens. 

Of  course  corrections  must  also  be  made  on  the  readings  of  the  transverse 
strain  due  to  the  effect  of  transverse  sensitivity  of  the  electric  strain 
gages.  The  transverse  sensitivity  factor  for  A-7  strain  gage  is  given  by 
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[3]  c 

Perry  and  Lissner  as  k  =  -0.01.  The  error,  which  is  defined  by  e  1  — - -  , 

where  c is  the  apparent  strain  and  c  is  the  true  strain,  is  3%  for  a  Poisson’s 

ratio  of  0.3  and  is  2.6%  for  a  Poisson's  ratio  of  0,33.  It  is  equivalent  to 

saying  that  e  =  0.971  ec  for  material  with  Poisson's  ratio  of  0.3  and  e  = 

0.984  c  for  material  with  Poisson's  ratio  of  0.33.  We  have  thus  seen  that 
c 

the  true  transverse  strain  is  2.9%  smaller  than  the  apparent  transverse 
strain  when  v  =  0.30,  and  is  1.6%  smaller  than  the  apparent  transverse  strain 
when  v  =  0.33.  This  correction  is  however  in  a  different  direction  from  the 
one  due  to  the  misalignment  of  the  strain  gages,  therefore  some  of  the  errors 
should  cancel  out. 

[4] 

Stang,  Greenspan  and  Newman  reported  in  as  early  as  1946  their  experi¬ 
mental  study  on  Poisson's  ratio  of  aluminum  alloys  24  ST  and  24  SRT,  chrome- 
molybdenum  steel  plate  and  structural  and  fully  killed  low-carbon  steel  plate. 
In  all  cases  these  authors  reported  the  increase  of  Poisson's  ratio  beyond 
its  initial  (elastic)  value  throughout  the  longitudinal  strain  range  tested. 
This  observation  is  not  in  agreement  with  the  present  findings. 

We  like  to  mention  that  the  above  discussed  misalignment  in  the  trans- 

[4] 

verse  strain  gages  also  existed  in  the  tests  by  Stang  et  al  .  In  addition 
to  this,  as  it  was  pointed  out  by  these  authors  themselves  in  their  report, 
there  were  large  discrepancies  in  the  values  of  Poisson's  ratio  obtained  by 
them  for  two  nominally  identical  specimens. 

Finally,  from  a  thermodynamic  point  of  view  the  hydrostatic  stress  versus 
volumetric  strain  curve  should  bend  over  toward  the  e^-axis^;  this  agrees 
with  our  own  observations. 
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VI .  Conclusions 

The  following  conclusions  can  be  drawn  from  the  present  experimental 
study: 

(1)  The  curve  of  hydrostatic  stress  versus  volumetric  strain  bends 
over  toward  the  volumetric  strain  axis  for  all  three  materials 
tested.  This  implies  that  the  commercially  pure  aluminum,  the 
electrolytic  tough  pitch  copper  and  the  low-carbon  steel  are  plas¬ 
tically  compressible,  (at  least  when  the  hydrostatic  stress  is  of 
a  tensile  character)  and  that  the  most  important  assumption  in  the 
classical  theory  of  plasticity  concerning  the  plastic  incompress¬ 
ibility  of  material  lacks  experimental  justification. 

(2)  Lateral  instability  occurs  weakly  in  the  case  of  aluminum  and 
occurs  strongly  in  the  case  of  low-carbon  steel.  Lateral  instabil¬ 
ity  always  follows  the  longitudinal  plastic  flow  for  aluminum, 
whereas  two  possibilities  arise  in  the  case  of  low-carbon  steel- 
lateral  instability  may  precede  or  follow  the  longitudinal  plastic 
flow  as  dictated  by  the  anisotropy  of  the  material.  No  phenomenon 
of  lateral  instability  has  been  observed  for  the  copper  specimens 
tested. 

(3)  The  lateral  instability  corresponds  to  a  decrease  in  volume  due  to 
a  slight  increase  in  load  for  the  case  of  aluminum,  and  corresponds 
to  a  sudden  decrease  in  volume  while  the  load  remains  constant  for 
the  case  of  low-carbon  steel.  Where  the  lateral  instability  occurs, 
Poisson's  ratio  has  a  large  increase. 

(4)  Poisson's  ratio  for  the  commercially  pure  aluminum  decreases  consid¬ 
erably  after  yield  has  occurred,  from  an  initial  value  of  0.30  to  a 


15 


value  as  low  as  0.06  corresponding  to  a  longitudinal  strain  of 
approximately  2%;  beyond  this  point  it  increases  continuously 
and  at  an  axial  strain  of  3%  it  reaches  a  value  of  0,15  which  re¬ 
mains  approximately  constant  up  to  =  9%. 

(5)  For  the  electrolytic  copper,  Poisson’s  ratio  decreases  continuously 
with  c  from  a  value  of  0.35  to  a  value  of  0,3  at  c  =  9%. 

X  X 

(6)  For  the  low-carbon  steel,  two  possibilities  prevail.  If  the  later¬ 
al  instability  occurs  before  the  longitudinal  plastic  flow  then 
Poisson’s  ratio  increases  suddenly  at  yielding  from  its  initial 
value  of  0.29  to  a  value  as  high  as  6.0.  It  then  decreases  grad¬ 
ually  and  reaches  a  value  of  0.3  at  -  3%,  Beyond  this  point, 
Poisson's  ratio  remains  approximately  constant  up  to  c =  9l.  If 

the  lateral  instability  occurs  after  the  longitudinal  plastic  flow, 

% 

then  Poisson’s  ratio  decreases  considerably  after  yield  has  occurred 
from  an  initial  value  of  0.29  to  a  value  as  low  as  0.02  correspond¬ 
ing  to  a  longitudinal  strain  of  approximately  2.5%,  beyond  this 
point,  it  increases  continuously  and  at  an  axial  strain  of  3%  it 
reaches  a  value  of  0,3  which  remains  approximately  constant  there- 
after. 

(7)  Loop3  similar  to  the  hysteresis  loops  are  present  in  curves  obtained 
by  plotting  -c^  against  and  also  in  curves  of  against 

for  all  three  materials  tested. 
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Abstract 

The  endochronic  theory  of  plasticity,  developed  previously  by  the  first 
author,  has  been  shown  to  give  an  accurate  prediction  of  metallic  response  in 
the  plastic  range.  The  mechanical  behavior  of  an  isotropic  material  is  given  in 
terms  of  two  "heredity"  functions  similar  to  those  encountered  in  the  theory  of 
linear  viscoelasticity.  In  the  present  paper  the  "tensile"  and  "transverse" 
heredity  functions  of  copper  and  aluminum  are  determined  experimentally  using 
the  above  theory.  The  properties  of  these  functions  are  also  discussed. 


1. 


1 .  Jn t roduct  ion 

The  recently  proposed  "endochronic  theory  of  plasticity"  has  been  shown 
to  give  accurate  analytical  predictions  to  a  variety  of  experimental  observa¬ 
tions  on  metals  such  as  copper  and  aluminum^.  This  theory  differs  from 
previous  theories  in  that  a  yield  surface  is  not  necessary  for  the  development 
of  the  constitutive  equations  and  that  the  stress  (under  isothermal  conditions) 
is  determined  by  the  previous  deformation  history  defined  on  a  time  scale  which 
is  independent  of  the  real  time,  but  is  itself  a  property  of  the  material  at 
hand. 

In  a  previous  paper^  it  was  shewn  that  diverse  phenomena  such  as  cross¬ 
hardening*,  loading-unloading  loops,  cyclic  hardening  and  effects  of  prestress, 
can  be  described  quantitatively  and  accurately  with  a  single  constitutive  equa¬ 
tion. 


Specifically,  this  equation  has  the  following  form  for  small  deformation 
and  initially  isotropic  materials: 
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*  This  term  means  a  change  in  the  uniaxial  stress-strain  relation  due  to  a 
torsional  prestrain  or  vice-versa. 


2. 


The  symbol  z  denotes  a  positive  monotonical ly  increasing  time  scale  wi th 
respect  to  a  time  measure  dr,  such  that 


'  pijk*  dE,j  dEki 
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where  is  a  material  tensor,  which  is  positive  definite  and  which,  for 

the  isotropic  materials  envisioned  in  Eq.  (l.l),  has  the  form 
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(1.6) 


where 


kj  and  k^  are  material  parameters,  such  that  kj  +  —  >  0,  k^  >  0. 


Materials  described  by  Eq.  (1.1)  are  strain  history  dependent  but  strain- 
rate  independent.  The  derivation  of  constitutive  equation  (1.1)  was  given  in 
deta i 1  in  Ref.  1 . 


It  may  be  seen  that  eq.  (l.l)  Is  uniquely  described  by  means  of  two 
material  "heredity  functions1*,  X(z)  and  p(z).  However,  eq.  (l.l)  may  be  • 
described  instead  in  terms  of  the  function  p(z)  and  K(z),  the  "shear  end  bulk 
heredity  functions",  respectively.  Specifically: 
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and  Sjj  and  Cjj  are  the  stress  and  strain  deviaters.  The  functions  A(z), 
p(z)  and  K(z)  are  related  as  follows: 


K(z)  »  A  (z)  +  j  w(z) 


(1.8) 


From  an  experimental  viewpoint  it  is  more  convenient  to  relate  *uw  and 
Cjj(z)  to  the  functions  E(z)  and  v(z),  which  we  shall  call  the  "tensile"  and 
"transverse"  heredity  functions.  This  relationship  is: 
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The  functions  E(z)  and  v(z)  are  iclated  to  the  previous  functions  through  their 


Laplace  transforms  as  follows: 
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where  a  bar  over  a  function  denotes  its  Laplace  transform  with  respect  to  the 
parameter  p. 


In  the  present  paper  we  present  forms  of  the  functions  E(z),  v(z)  that 
have  been  determined  experimentally  in  the  case  of  copper  and  aluminum;  these 
will  be  discussed  in  later  sections.  The  experimental  measurements  necessary 
to  determine  the  above  functions  were  also  used  to  evaluate  critically  the 
implications  of  the  assumption  of  plastic  incompressibility,  normally  made  in 
the  classical  theory  of  plasticity.  These  experiments  are  reported  in  detail 
In  Ref.  2. 


As  Is  well  known,  this  assumption  leads  to  the  result  that  Poisson's  ratio 
tends  to  y  as  the  plastic  strain*  increases.  Our  experimental  results  do  not 
support  this  conclusion.  In  the  classical  theory  of  plasticity,  plastic  in¬ 
compressibility  has  also  been  interpreted  to  mean  that  the  hydrostatic  stress 

*  The  term  "plastic  strain"  is  used  in  the  content  of  the  classical  theory 
of  plasticity.  The  endochronlc  theory  does  not  recognize  the  dichotomy 
of  strain  into  elastic  and  plastic  parts. 


i 


k. 


is  a  linear  function  of  the  hydrostatic  strain;  our  experimental  data  do  not 
support  this  either. 

Of  course  in  its  broader  aspect,  plastic  incompressibility  may  imply  a 
reversible,  though  non-linear,  volumetric  response.  However,  the  observed 
volumetric  response  was  decidedly  irreversible.  Consequently  our  experimental 
measurements  do  not  support  the  assumption  of  plastic  incompressibility,  at 
least  when  the  hydrostatic  stress  is  of  tensile  character.  (A  non-linear 
irreversible  volumetric  response  in  simple  compression  has  also  been  reported 
in  the  literature  by  Bridgman^.) 

It  is  significant  that  the  endochronic  theory  can  predict  plastic  com¬ 
pressibility,  exactly,  provided  that  the  heredity  function  K(z)  is  chosen 
appropriately,  or  is  determined  experimentally. 


In  Ref.  1,  the  time  scale  z  was  related  to  the  intrinsic  time  C  by  the  equation 

2  *  j  log  (1  +  B?)  (2.3) 

which  proved  very  satisfactory  in  providing  agreement  between  the  theory  and  a 
variety  of  experiments  as  discussed  in  the  Introduction.  For  this  reason  eq. 
(2.3)  will  be  retained  In  this  paper. 


5. 


As  a  result  of  eq.'s  (1.5)  and  (1.6) 

d?2  =  k,2 (dex  +  2ey)2  +  k2(de2  +  2dey)  (2.0 

Because  kj  and  k^  are  not  known  and  to  illustrate  the  methodologv  f  the 
theory,  we  have  taken  k,  =0.  In  this  event 

1/2 

dex  (2.5) 

If,  in  addition,  an  experimental  relation  between  a and  e^  as  well  a 

between  e  and  e  were  known,  then  from  eq.'s  (2.3)  and  (2.5)  e  would  be  a 
y  x  x 

known  function  of  z,  in  which  case  o  and  e  would  also  be  known  functions 

’  x  y 

of  z.  >• 

Equations  (2.1)  and  (2.2),  which  are  Volterra  integral  equations,  would 
then  be  solved  to  yield  the  functional  form  of  the  Kernels  (heredity  functions) 
E(z)  and  v(z). 


In  the  particular  use  where  is  a  linear  function  of  z  (experiments 
show  that  this  is  a  very  good  approximation  as  will  be  shown  later)  eq.'s  (2.1) 
and  (2.2)  can  be  solved  exactly  to  yield  the  relations: 
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(2.8a) 

In  Ref.  2  we  presented  experimentally  obtained  uniaxial  stress-strain 

* 

curves  for  a  large  number  of  copper  and  aluminum  specimens  .  We  also  present 

ed  experimental  curves  showing  the  relation  between  e  and  r  for  these  same 

y  ^ 


Electrolytic  tough  pitch  copper  and  commercially  pure  aluminum. 
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specimens.  In  Fig.'s  !-*»  we  show  such  typical  curves 

Now,  using  the  experimental  curves  of  Fig.'s  3  and  4,  eq.  (2.5)  was 
determined  as  a  function  of  £  for  copper  and  aluminum  and  the  curves  of  these 
functions  arc  shown  in  Fig.  5»  with  remaining  undetermined.  Similar  curves 
were  obtained  for  as  a  function  of  £  and  appear  in  Fig.  6. 

One  can  see  in  Fig.  5  that  the  relation  between  and  £  is  linear  to  a 
very  good  approximation  for  both  copper  and  aluminum.  Therefore,  eq.'s  (2.6) 
and  (2.7)  may  be  considered  as  approximate  solutions  of  the  Vol terra  integral 
equations  (2.1)  and  (2.2)  and  may  be  used  to  determine  the  form  of  the  hereditv 
functions  E(z)  and  v(z).  This  was  done  as  follows:  Eq.'s  (2.6)  and  (2.7)  may 
be  written  in  the  alternative  form: 

E(z)  »  e6z  -  HBo(z)  (2.9) 

*v(z)  *  57^  -  hBe  (z)  (2.10) 

°Cx  y 

where  now  h  Is  taken  for  consistency  to  be  the  mean  value  of  the  slope  of  the 
approximate  linear  -elation  between  and  5.  Using  eq.'s  (2.9)  and  (2.10) 
and  Fig.'s  (l-M»  E(z)  and  v(z)  wore  calculated  for  copper  and  aluminum;  these 
appear  in  Fig.'s  7*8  and  9*10  for  various  values  of  S. 

Discussion  of  the  heredity  functions 

(a)  Copper  Two  characteristics  of  the  tensile  heredity  function  E(z) 
merit  discussion.  Ir,  the  region  of  small  z,  E(z)  exhibits  a  point  of  inflexion 
this  Implies  that  the  second  derivative  of  E(z)  changes  sign  at  that  point. 

** 


For  copper  specimen  k  and  aluminum  specimen  k  of  Ref.  3. 
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This  behavior  is  not  predicted  by  the  current  form  of  the  theory,  according  to 
which 

n 

E  (z)  «  I  E  e‘V  (2.11) 

r=1  r 

<j2  £ 

where  Er  and  ar  are  positive  constants.  Evidently  ,  as  calculated  fron) 

the  above  equation,  is  always  positive.  However,  except  for  this  small  region 
in  the  vicinity  of  z= 0,  the  experimental  curve  assumes  a  shape  not  unlike  the 
one  predicted  by  the  eq.  (2.11).  On  the  other  hand  the  form  of  the  function 
E(z)  does  not  appear  to  be  particularly  sensitive  to  changes  in  8,  particularly 
at  moderately  small  values  of  z. 

The  transverse  heredity  function  v(z)  decreases  slowly  as  z  increases  but 
remains  virtually  constant  (over  the  range  of  2  shown)  if  8  Is  small.  There¬ 
fore,  as  far  as  copper  is  concerned  the  assumption  of  constant  Poisson's  ratio 
Is  true,  to  a  very  good  approximation. 

(b)  Aluminum:  The  tensile  heredity  function  for  aluminum  is  very  similar 
to  that  of  copper.  The  transverse  heredity  function,  however,  exhibits  a  strik¬ 
ingly  different  behavior.  At  moderate  strains  it  decreases  from  Its  initial 
value  quite  rapidly,  until  it  becomes  virtually  zero;  then,  at  higher  strains 
It  increases  abruptly  to  a  value  which  is  three  times  the  Initial  value.  Then 
it  decreases  again  and  at  still  higher  strains,  it  tends  tu  a  constant  value, 
which  is  nearly  equal  to  the  initial  value. 

No  physical  explanation  can  be  offered  at  this  time  for  this  behavior, 
which  is  certainly  not  due  to  necking  of  the  specimen;  quite  definitely,  neck¬ 
ing  did  not  occur  in  thr*  range  of  deformation  indr"'ted  in  Fig.  8. 


3 .  Mat  e  rial  compress  ib  i  1  i  ty  u  n  d  or  hydros,  t  at  ic  t  c  n  sj  on 

In  this  section  we  shall  discuss  the  experimental  measurements  which  led 
to  the  conclusion  that  the  assumption  cf  plastic  incompressibility  (normally 
made  in  the  classical  theories  of  plasticity)  is  not  tenable. 

The  experiments,  simple  in  nature,  consisted  in  monitoring  the  axial 

stress  as  well  as  the  axial  strain  and  one  transverse  strain  of  a  thin  flat 

* 

bar  In  uniform  uniaxial  tension  .  This  simple  test  has  enjoyed  a  great  deal 

(i,  i 

of  popularity  but  we  could  find  only  one  Instance  '  where  "transverse  strain" 
measurements  were  made  in  the  plastic  range.  This  is  surprising,  since  this 
test  affords  a  very  critical  evaluation  of  the  assumption  of  plastic  incom- 
pressfbll Ity. 

The  classical  plasticity  theory  yields  the  following  result,  assuming 
elastic  comprcsstbil i ty  and  a  Von  Mises  yield  surface: 

rfcY  1  I 

-  -  V  S  <3  ’> 

where  e  and  e  are  the  longitudinal  and  transverse  strains  respectively, 
and  are  the  clastic  Young's  modulus  and  Poisson's  ratio  respectively, 
and  E^  Is  the  tangent  modulus  of  the  tensile  ,;trcss-strain  curve. 

£q.  (3-D  can  be  Integrated  numcr ical ly  to  yield  a  relation  between 
and  c^.  This  has  been  done  for  aluminum  and  copper. 

The  aluminum  specimens  were  sheared  from  a  1100-0  aluminum  sheet  and 
annealed  at  600°f  for  two  hours  and  then  oven  cooled  to  room  temperature. 

The  copper  specimens  were  cut  out  from  an  electrolytic  tough  pitch  copper 
110(99-9  4  X  cu)  bus  bar  and  annealed  at  75u°f  and  then  even  cooled  to  rc<x' 
temperature. 


Detailed  description  of  the  apparatus  and  the  experimental  p  occdn-e  Is  to 
be  found  in  Ref.  2. 


he  results  are  shown  in  Fig.'s  3  and  b.  It  may  be  seen  that  the 
classical  plasticity  theory  gives  a  poor  prediction  of  transverse  strain  re¬ 
sponse,  particularly  for  aluminum. 

Cn  the  other  hand  the  assumption  of  constant  Poisson's  ratio  is  quite 
good  over  a  wide  range  of  axial  strain  for  aluminum  and  is  almost  exactly 
true  in  the  case  of  copper.  For  this  reason  the  assumption  of  constant 
Poisson's  ratio  was  adopted  in  the  analysis. 

Also  shown  in  Fig.'s  11  and  12  are  plots  of  cr  ^  versus  e  ^  These  plots 
show  very  vividly  that  hydrostatic  response  is  not  elastic,  at  least  when  the 
hydrostatic  stress  is  tensile.  If  one  extrapolates  from  one's  experience  in 
elasticity,  the  effect  of  this  assumption  is  likely  to  be  large  in  the  case 
of  kinematically  constrained  configurations. 
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